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The determination of the molecular weights of large molecules by measuring the turbidity of 
the solution and the change in index of refraction with concentration is discussed. The appa- 
ratus, its calibration, and the technique used are treated togethér with a comparison between 
molecular weights measured by this method and other methods. The effect of polymolecularity 


is also presented. 





INTRODUCTION 


ECENTLY Debye! has pointed out that the 
measurement of light scattering in solutions 
may be used to determine the molecular weights 
of high polymers. Since that time we have been 
interested in examining this method and some 
data of preliminary experiments have already 
been published.? In the continuation of this in- 
vestigation a few anomalies have been noted 
which are as yet incompletely understood. How- 
ever, in view of the apparent success of this 
method, we report here in more detail concerning 
the experimental methods used and some further 
results obtained. Some limitations and pitfalls as 
well as an extension of the method to the charac- 
terization of molecular weight distributions are 
discussed. 
Before proceeding we review briefly the theo- 
retical background.** It is common knowledge 
1P. Debye, J. Apps Phys. 15, 456 (1944). 
2 P.M. Doty, B. H. Zimm, and H. Mark, J. Chem. Phys. 
12, 144 (1944). 
3S. Bhagavantam, The Scattering of Light and Raman 
Effect (Anhandra University, Waltair, India, 1940); 
American Edition (Chemical Publishing Company, Brook- 
lyn, New York, 1942). 


4J. Cabannes, La Diffusion Moléculaire de la Lumiére 
(Les Presses Universitaires de France, Paris, 1929). 


that some light is always scattered from a beam 
of light passing through a transparent fluid 
(Tyndall effect). This, however, would not occur 
in a medium which was perfectly uniform. Real 
solutions are never perfectly uniform because of 
minute thermal fluctuations in density and con- 
centration. These inhomogeneities are the source 
of the scattered light. The extent of the fluctua- 
tions is determined by the free energy change 
(work) necessary to produce them. On the basis 
of these considerations, von Smoluchowski® and 
Einstein® derived the following expression for the 
turbidity of a colorless binary mixture. 
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to, the turbidity of the solution, corresponds 
to the extinction coefficient in Lambert’s law, 


5M. von Smoluchowski, Ann. d. Physik 25, 205 (1908). 
6 A. Einstein, Ann. d. Physik 33, 1275 (1910). 
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i.e., it is the fractional decrease in the intensity 
of incident light upon transversing one centi- 
meter of solution, the decrease occurring only 
from scattering—not absorption. The other 
symbols are: \= wave-length in vacuum (cm~) ; 
No=Avogadro’s number ; R= molar gas constant ; 
T =temperature (°K) ; 8=isothermal compressi- 
bility; p=density; »=refractive index of solu- 
tion; M,=molecular weight of solvent; f1=fu- 
gacity of solvent; and c=concentration of solute 
in grams per cm? of solution. 

While Eq. (1) has been used successfully in this 
and previous investigations in this laboratory, 
and is essentially the one used by Debye,} it 
should be remarked that there is not unanimous 
agreement on the value of the constant factor 
before the braces in Eq. (1). Several authors have 
proposed the introduction of a factor [(m?+-2)/3 ? 
in this constant. The arguments are well sum- 
marized by Cabannes.‘ Since previous experi- 
ments‘ designed to decide the question (by 
measurement of the light scattered by one-com- 
ponent. liquids) have apparently been incon- 
clusive, we have not yet attempted to decide the 
question but have used the formula as originally 
proposed. 

The first term in Eq. (1) represents the 
turbidity arising from fluctuations of density. 
This part of 79 corresponds essentially to the 
turbidity of the pure solvent (assuming 8 to be 
constant). Thus if the measured turbidity of the 
solvent is subtracted from 7» the remainder, 1, 
- represents the turbidity caused by concentration 
fluctuations, a quantity directly related to the 
solute molecular weight. 

Substituting in (1) the thermodynamic relation 


—@ In f1/ac=(Vi/RT)(dx/dc), (2) 


TABLE I. Comparison of molecular weight measurement. 








Molecular weight 


Light 
scattering 


178,000* 
182,000 
107,000 
190,000 
445,000 
123,000 
76,000 


Osmotic 
pressure 


172,000 
198,000 
91,000 


Material 





Polystyrene, No. 
Polystyrene, No. 
Polystyrene, No. 
Polystyrene, No. 
Polystyrene, No. — 
Cellulose Acetate, No. 1 (125,000) 
Cellulose Acetate, No. 3 75,000 








* This value was checked independently by two other laboratories. 
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there is obtained (2) 
322°cn?(dn/dc)? 
T= ’ 
3\4No(1/RT)(dx/dc) 





(3) 


where = is the osmotic pressure and V the partial 
molal volume of the solvent has been replaced by 
its approximate equivalent M;/p. The osmotic 
pressure of high-polymer solution of low concen- 


tration is given by 
a=(RT/M:2)c+Be?+:--, (4) 


which, when substituted into (3), gives upon 


rearrangement 
c 3r4No fi 2B 
-= —+—et «| 
+ 32e%n%(dn/ac)LM2 RT 
1p 1 2B 


dl —+—et--| (5) 
M,. RT 





Thus it appears that a plot of Hc/r versus c 
should give at low concentration a straight line 
whose intercept is the reciprocal of the molecular 
weight of the solute. Furthermore, the slope of 
this plot is twice the slope obtained from osmotic 
pressure measurements when 2/cRT is plotted 
against c. In this way the measurement of the 
two optical quantities, (1) turbidity as a function 
of concentration and (2) the change of index of 
refraction with concentration, gives information 
similar to that obtained from direct osmotic 
pressure measurements. 

However, it is important to note that whereas 
osmotic pressure measurements yield a number 
average molecular weight, turbidity measure- 
ments yield a weight average.* In the final section 
a more exact discussion of this point occurs, but 
the conclusion can be seen intuitively from a 
molecular picture of light scattering. In a very 
dilute solution the intensity of the scattered 
light is equal to the sum of the intensities of the 
light scattered by the individual molecules.’ The 
intensity of the light scattered by a single mole- 
cule of a given species is proportional to the 
square of the molecular weight. Since the number 
of molecules is equal to the total weight of that 

* The number average molecular weight (2M; Ni/2Ni) 
counts each molecule with equal weight, whereas the weight 


average (2M;*N;/2M;:N;) weights each molecule ac- 
cording to its mass. 
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Fic. 1. Plot of turbidity data for polystyrene fractions in 
methyl ethyl ketone. 


species divided by the weight of the individual 
molecule, the total intensity of light scattered 
from one species is therefore proportional to the 
weight of that species multiplied by its molecular 
weight. Averaging over all molecular weight 
species yields the weight average molecular 
weight. This conclusion has been experimentally 
demonstrated.’ 


EXPERIMENTAL DETAILS 


Numerous methods may be adapted to measure 
the turbidity of high-polymer solutions. One 
could hope to measure 7 directly by measuring 
the transmission of the solutions. This is usually 
unsatisfactory, however, because of the small 
values of 7 (10-? to 10-* cm~") and the frequent 
presence of absorbing impurities. On the other 
hand, direct measurement of the intensity of the 
scattered light conveniently avoids these diff- 
culties and readily lends itself to visual or photo- 
electric observation. Moreover, it has been our 
observation that transmission measurements are 
subject to a number of uncertainties, the absence 
of which can be assured only by measurements 
over a rather large spectral range. True absorp- 
tion, low angle scattering, small amounts of 
impurities, solute size comparable with wave- 
length, and other effects still not understood, may 


7B. H. Zimm and P. M. Doty, J. Chem. Phys. 12, 203 
(1944), 


{3 


















































5 





0 4 CONC. 8 1.2 


Fic. 2. Plot of turbidity data for polystyrene fractions in 
methyl ethyl ketone. 


prevent the turbidity as measured by transmis- 
sion from being the proper value to use in Eq. (5). 

Although photoelectric methods are always 
capable of greater precision, the visual method 
has certain inherent advantages, such as small 
aperture, direct view of scattering solution, exact 
linearity, and freedom from electrical difficulties. 
We have constructed a turbidimeter attach- 
ment** for use with the Zeiss Pulfrich photome- 
ter. This turbidimeter was patterned after the 
Zeiss turbidimeter, except that the light scattered 
is observed at an angle of 90° to the incident 
beam, instead of 135°. Light from a mercury bulb 
(AH-4 100 watt G.E.) is rendered parallel by lens 
and passes through the solution in a square cell. 
A small, constant fraction of the incident light is 
reflected from a glass plate (mounted in the 
beam) to an opalescent glass which serves as a 
standard comparison source for the photometer. 
The dial reading corresponding to matched fields 
in the photometer is directly proportional to the 
total turbidity. 

It is necessary to calibrate the instrument; that 
is, to determine the proportionality constant 
between dial reading and absolute turbidity. One 
method of calibration is based on measuring the 
ratio of the intensities of the incident and 


** General Gauge and Tool Company, 30 Irving Place, 
New York. 
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scattered light, with the help of a white (mag- 
nesium carbonate) surface making a 45° angle 
with the incident beam in place of the solutions. 
If the solution and the magnesium carbonate are 
measured successively, the ratio of the photometer 
readings gives the quantity 


(ir?/I)(ha/sin 45°), 


where 7 is the intensity of the light scattered from 
one cubic centimeter of solution at distance 7, I is 
the intensity of the incident beam, and h is the 
width (in the direction of observation) of the 
incident beam of homogeneous density. If the 
scattering of the solutions obeys Rayleigh’s law 
for its angular distribution, i.e., if it varies as 
(1++-cos? 6) where @ is the angle between incident 
beam and direction of observation, integration 
gives 


ro= (16/3) (ir?/I). 


Correction must be made for the reflectivity of 
magnesium carbonate (97.5 percent) and for re- 
flection at the glass-air interfaces. Thus the 
calibration consists of the determination of the 
ratio of the photometer readings with solution 
and block in the cell successively together with 
the measurement of the width of the incident 
beam in the cell. 

Another means of calibration is based on the 
use of standard opalescent glass prisms furnished 
by Carl Zeiss and Company. A polymer solution 
was measured in a Zeiss turbidimeter using two 
of these standard prisms. The value of the 
turbidity so obtained checked satisfactorily the 
turbidity found for the samé solution by the 
previous method. 

The measurement of the transmission of 
a strongly scattering colorless solution offers 
another method of determining absolute tur- 
bidity. Several preliminary attempts yielded 
inconsistent turbidity values that were also too 
large, probably indicating absorption. The in- 
vestigation of this method is continuing. 

The sensitivity of our instrument is such that 
To values as small as 10~* cm~! may be measured 
with a reproducibility of 2 percent. 

’The other optical measurement necessary is 
the evaluation of dn/dc. Our measurements have 
confirmed that » varies linearly with c at low 
concentrations. Consequently one measurement 














ZIMM, AND H. MARK 





of the refractive index difference between the 
solvent and solution determines dn/dc. Since 
usually the magnitude of the difference for a 1 
percent solution is approximately 0.001, a differ- 
ential method is desirable. We have used a 
Pulfrich refractometer with a divided cell for this 
purpose, thus permitting a precision of about 1 
percent in this measurement. 

Depolarization (ratio of the intensities of the 
horizontal to vertical components in the scattered 
light) can easily be measured by the use of a 
Polaroid or nicol inserted in the solution side of 
the photometer. Strictly speaking, Eqs. (1) and 
(5) should contain a factor (Cabannes), depend- 
ing on the depolarization, which arises from 
fluctuations in orientation of anisotropic mole- 
cules.'* 4 In the work reported here the de- 
polarization was so small that this factor was 
essentially unity. 

All solutions and solvents were filtered through 
fine or medium sintered glass filters. A deep glass 
cell (33X10 cm) has been used in later work. 
This permits one to begin measurements with 
25 cc of (say) a 1 percent solution. Measured 
amounts of filtered solvent may be added and 
stirred; thus all the measurements at different 
concentrations are made in one continuous 
operation. 


RESULTS 


In this section typical examples of the use of 
this method on some fractionated high polymers 
are presented together with available osmotic 
pressure data. The results of the turbidity meas- 
urements and osmotic pressure measurements are 
shown in Table I for the fractions studied. 

The most carefully measured fraction (No. 1) 
was one prepared and measured osmotically by 
Hanson.* The other two polystyrene fractions 


(Nos. 2 and 3) were prepared and their osmotic 


pressure determined in this laboratory several 
years ago. In Figs. 1 and 2 the turbidity data 
are plotted in such a manner that the intercept is 
the reciprocal of the molecular weight. The data 
for two other fractions for which no other mo- 
lecular weight data are available are also shown 
in Fig. 1. 

The need for extrapolation of H(c/r) to infinite 


8 J. Hanson, Thesis, Polytechnic Institute of Brooklyn, 
1944, 
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dilution is clearly shown in the figures. For ex- 
ample, if the molecular weight of sample No. 5 
were calculated from the turbidity of its 1 


percent solution alone, the value would have been’ 


110,000, instead of the actual value of 445,000. 
The higher the molecular weight, the larger the 
error which results from use of data without 
extrapolation. 

The cellulose acetate fractions were prepared 
by Badgley.® The molecular weight in parenthesis 
is from viscosity measurements. However, the 
viscosity-molecular weight equation for this 
series of fractions has been so well established® 
that the value has the validity of an osmotic 
pressure measurement. 

It is interesting to examine the relation be- 
tween the slope of the H(c/r) versus c plot and 
the w/c versus c plot. As stated previously, if 
osmotic pressure is represented by 


a=(RT/M:2)c+Bce?, (4) 


the slope of the turbidity graphs should give the 
value of 2B/RT. 
If one uses Huggins’ equation’® 


« RIT RT 
mf 
Cc M2 Vid,? 





(3 —p)¢, (6) 


where V; is the molar volume of solvent, dz is 
the density of solute, 4, the empirical constant 
of the solvent-solute system, may be used to 
characterize the slope, B, of the r/c versus c plots. 


That is: 
A(r/c) RT , 


+71 \ 
Ac V,d_? 





2—H1), 


TABLE II. Comparison of yu values. 











p» (Osmotic » (Light 
pressure) scattering) 
Polystyrene-toluene 0.43-0.45%" 0.43-0.44¢ 
0.44 
Polystyrene—MEK — 0.47-0.48* 
Cellulose acetate-acetone 0.42-0.43¢ 0.44* 








* This research. 
ons} Alfrey, A. Bartovics, and H. Mark, J. Am. Chem. Soc. 65, 2319 


(1 


> Reference 8. 
¢ Reference 2. 
4G. V. Schultz, Zeits. f. physik. Chemie A176, 317 (1936). 
¢ Reference 9. 


® W. Badgley, Thesis, Polytechnic Institute of Brooklyn, 
10M. L. Huggins, Ann. N. Y. Acad. Sci. 41, 1 (1942). 
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and 
ACH(c/7)]_ 2 


— \ 
Ac Vid 2” 





2 — M1). (7) 


In other words, it is convenient to compare the 
slopes of the lines obtained from plots of osmotic 
pressure and turbidity data by calculating the 
ui, from Eq. (7). In Table II the yu; values of 
the polystyrene-toluene system (previously re- 
ported) and cellulose acetate-acetone system 
obtained from turbidity data are compared with 
ui calculated from osmotic pressure measure- 
ments on the same system. 


DISCUSSION 


When the particle size of the solute becomes 
comparable to the wave-length of light, interfer- 
ence between waves from different parts of the 
same particle occurs and the simple theory 
outlined above breaks down. The total scattered 
intensity then falls below that given by Eq. (5). 
At the same time, instead of the angular distri- 
bution of the scattered light being symmetrical 
about 90° (as is the case with small particles), 
interference in the case of large particles causes 
more light to be scattered in the forward than in 
the backward direction. Measurement of the 
symmetry of the scattered intensity therefore 
furnishes a test of the applicability of Eq. (5). 

For spherical particles, deviations from Eq. (5) 
become noticeable when the particle diameter 
exceeds about one-tenth the wave-length of light. 
Polymer chain lengths frequently are greater 
than this by a large margin so that it might be 
expected that interference would often occur. 
However, most high-polymer molecules are be- 
lieved to be loosely convoluted in solution so that 
the longest dimension of the molecule is on the 
average much less than the chain length. It is 
clear that this point is critical, for the average 
extension of a randomly coiled polymer molecule 
of moderate weight will not exceed one-tenth the 
wave-length of light, whereas the common inter- 
pretation of viscosity measurements, sedimenta- 
tion constants, and flow birefringence indicate 
that most polymer molecules are more extended 
in solution and hence should give rise to a 
dissymmetry of the angular intensity of the 
scattered light. 

A simple instrument has been constructed to 
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measure the dissymmetry. A parallel beam of 
monochromatic light is passed through the solu- 
tion in a five sided cell. The light scattered at 45° 
on each side of the normal plane is reflected by 
the mirrors to the photometer apertures. The 
relative intensities are then compared by means 
of the photometer. 

Measurements of this type on the solutions of 
polystyrene in methyl ethyl ketone showed 
practically no dissymmetry, thus permitting the 
direct application of Eq. (5). Incidentally the 
higher molecular weight polystyrene fraction 
No. 5, dissolved in toluene did exhibit a con- 
siderable dissymmetry indicating that the mole- 
cule is more extended in the better solvent. On 
the other hand, all of the cellulose acetates re- 
ported here gave noticeable dissymmetries in 
acetone solution. Corrections based on these 
measured dissymmetries range upward to 20 
percent and have been applied to the results 
listed in Table II. A detailed discussion of these 
measurements and corrections will be given later 
when this study of cellulose acetate fractions is 
finished. 

In most cases, solutions in the concentration 
range between 1.0 gram and 10.0 grams of 
polymer. per liter have been used both for 
turbidity and refractive index measurements. 
The polymer specimens were the result of re- 
peated fractionation of a commercial material or 
of a sample polymerized under certain controlled 
conditions. Two types of impurities, which can 
never be completely eliminated, require special 
attention and care: 

A. Low molecular weight: impurities of the 
solvent or the polymer. All solvents were carefully 
distilled until their characteristics, such as den- 
sity, boiling point, and refractive index, did not 
undergo any changes upon further purification. 
Some solvents, such as dioxane, formic acid, etc. 
retain water very tenaciously and have to be kept 
carefully dry, because even small amounts of 
water affect the refractive index noticeably. The 
polymers should be dissolved and reprecipitated 
before they can be brought into the solution 
which is finally used for the measurements. It 
seems advisable to prepare a 2-3 percent solution 
of the original polymer in the same solvent in 
which the final measurements are contemplated, 
and to precipitate it with a volatile precipitant. 
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The precipitate is filtered off and thoroughly 
washed with the precipitant. It is then slowly 
dried at low temperatures in vacuum. Fast drying 
under more vigorous conditions, such as elevated 
temperature, leads in certain cases to hornifica- 
tion of the samples, which makes subsequent 
dissolving difficult and eventually even leads to 
solutions containing minute gel particles, which 
may lead to excessive scattering. This process of 
purification through precipitation should be re- 
peated one or more times; some polymers, e.g., 
polystyrene, are relatively easily separated from 
traces of monomer, catalysts, or stabilizers, others 
e.g., Nylon or rubbery copolymers, retain im- 
purities rather tenaciously. It does not seem easy 
to prepare any strict prescription at present, and 
it will be necessary to collect further experience 
until certain standard procedures can be rec- 
ommended. 

B. Coarse impurities of the solvent and 
polymer. All solvents and solutions have to be 
filtered repeatedly until they appear optically 
clean and until further filtration does not de- 
crease the turbidity any more. All solvents con- 
tain a certain amount of dust of mostly inorganic 
origin, which can be satisfactorily reduced to an 
insignificant amount by repeated filtration 
through glass fritters. All polymers contain coarse 
impurities of various types. Dust of inorganic 
origin, having large particle sizes (several mi- 
crons), is comparatively easily removed by pre- 
cipitation, washing, redissolution, and final filtra- 
tion through glass fritters. Haze, of organic 
origin, is frequently owing to chemical and 
structural heterogeneity of the polymer, such as 
too highly (or too lowly) acetylated fractions in 
cellulose acetate or branched (even cross-linked) 
constituents in vinyl polymers and copolymers. 
This type of haze, which is chemically closely 
related to the investigated material and differs 
from it only slightly in its solubility charac- 
teristics, sometimes offers very serious obstacles 
and eventually can prevent the application of the 
light scattering method altogether. It seems that 
a comparatively successful way to reduce this 
haze is to precipitate the polymer three or four 
times, using each time another solvent and 
eventually also more than one precipitant. It can 
be directly observed in such cases that the 
turbidity of the solution decreases visibly after 














each precipitation, but even a solution which 
appears completely transparent to the eye may 
still contain enough haze particles to give wrong 
results. These very large complexes affect the 
light scattering in about the same way as they 
affect specific viscosity, because of the fact that 
they contribute to the total effect according to 
their weight. Their large size (1000A or more), on 
the other hand, makes it possible to obtain inde- 
pendent experimental evidence about their pres- 
ence by measuring the dependence of the turbidity 
upon the angle between incident and scattered 
light in the apparatus described previously. 

In choosing a solvent for this type of work, one 
must be mindful of the great importance of the 
difference in index of refraction between solvent 
and solute. This difference should be greater than 
0.0005 refractive index unit between a 1 percent 
solution and the solvent. In fact, it is just this 
difference which determines the lowest molecular 
weight which can be measured with a given 
accuracy, for, obviously, the greater the refractive 
index difference, the lower the molecular weight 
solute that will give the same scattering. The 
lowest molecular weights we have measured were 
near 10,000 where 0n/dc=0.07 for the system 
investigated. To measure molecular weights lower 
than this dn/dc would have to be increased by 
use of a different solvent or less accuracy would 
have to be expected. Thus at present the method 
is limited to molecular weights above 10,000. 

Colored solutions offer some difficulty; how- 
ever, by using light of a wave-length that is 
outside the absorption band of the solution no 
trouble is encountered. 

At the end of this section it may be reiterated 
that the few recommendations made are far from 
being a well worked-out technique. They only 
represent the limited experience collected during 
a short period and will have to be considerably 
improved and expanded in the future. 


MOLECULAR WEIGHT DISTRIBUTION 


It was shown in a previous letter by Zimm and 
Doty’ that Eq. (5) can be generalized to include 
the case of a solute containing molecules of 
different degrees of polymerization. By the same 
methods used to derive Eq. (1), the turbidity of a 
solution of several different solutes in the same 
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solvent can be shown to be 
322° M yn? (dn/dc;)c; 


= ’ (8) 
3Nod4p i (—d@ In f;/dc;) 








where c; is the concentration (g/cc) of the ith 


‘ species of solute. 


If all the solute species are homologous poly- 
mers of the same monomer, differing only in 
molecular weight, two important simplifications 
can be made. In the first place, one might 
reasonably expect dn/dc to be independent of the 
molecular weight so that it can be factored out of 
the sum. In the second place, it is generally known 
that the coefficient of the second term in the 
expansion of the osmotic pressure in terms of 
concentration (Eq. (1)) is independent of mo- 
lecular weight in a homologous polymeric series. 
We might expect, therefore, that this coefficient 
would be constant in a solution of mixed homologs 
as in the following equation 


RT 
r=> —e+Be?+Ce+::-, (9) 
M; 
where B is independent of molecular weight. 

Introduction of these assumptions leads’ to the 
following expansion for H(c/r) 


0 )~atel'*(r) ).f 
tele) ). | 
“erll(@) ).! 

Gr) ).(Gr) Jew 


where M is the weight average molecular weight 
and ((AM/M)?),, and ((AM/M)*) are the second 
and third moments of the molecular weight 
distribution curves about the weight average 
molecular weight. This reduces to Eq. (5) when 
the moments are all zero, i.e., when the solute 
consists of one molecular weight species only. 
It can immediately be seen that the plot of 
H(c/r) for solutions of mixed molecular weight 
may have appreciable (downward) curvature. 
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Fic. 3. Turbidity measurements on cellulose acetate 
in acetone. 


Experimental data ona solution of unfractionated 
cellulose acetate shown in Fig. 3 illustrate this 
predicted curvature. We see here an observable 
effect of non-homogeneity of molecular weight 
which has no counterpart in the case of osmotic 
pressure. This difference arises from the same 
reasons that cause osmotic pressure to be related 
to number average and turbidity to weight 
average molecular weight. 

One way of utilizing this difference would be 
to plot H(c/r) versus c for the unfractionated 
polymer and draw the best curve through the 
data. Then either from turbidity measurements 
on a fraction or from osmotic pressure measure- 
ments on the unfractionated polymer one may 
obtain the slope of the H(c/r) versus c plot for a 
fraction. The difference between this slope and 
the initial slope of H(c/r) versus c for the unfrac- 


tionated polymer will give the value of the second 
moment, 


((AM/M)*)w, 


which characterizes the “width” of the distri- 
bution curve or the dispersion of molecular 
weights about the average value. By ordinary 
curve fitting methods the coefficient of the c? in 
Eq. (10) may be determined, from which the 
third moment, 


((AM/M)*)w, 


may be evaluated using the previously found 
value for the second moment. This third moment 
characterizes the asymmetry of the molecular 
weight distribution curve. The determination of 
higher moments would lead to more specific 
information about the distribution curve but at 
present this seems impracticable both experi- 
mentally and theoretically, for extremely accu- 
rate measurements would be required and the 
compliance of Eq. (10) with the system in- 
vestigated would have to be rigorously proven. A 
study of the application of this general method is 
being continued. 
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The absorption spectrum of aniline vapor has been photographed in the first and second order 
of a 3-m grating spectrograph. The strongest bands appeared as doublets and on the basis of 
these, a tentative analysis was made and compared with the available Raman data. 





INTRODUCTION 


HE analysis of the near ultraviolet absorp- 
tion spectrum of benzene! has shown that 
only the transitions involving non-totally sym- 
metric vibrations of the e,+ type are allowed. It 
has been further shown? that a modification of the 
electronic structure by substitution tends to re- 
move the restriction so as to allow the 0—0 
transition and transitions not involving the e,* 
type of vibration. The modification of the 
electronic structure of the benzene nucleus occurs 
by virtue of an inductive effect which is depend- 
ent on the electronegativity of the substituent 
and a resonance effect which requires an unshared 
pair of electrons on the substituent. 

Sponer and Wollman* have shown the presence 
of two spectra, in monochlorobenzene one in- 
volving an ¢,+ vibration and other involving 
symmetric vibrations only, in addition to the 
0—0 transition. It was expected that the aniline 
spectrum would show the same features with the 
symmetrical frequencies appearing more strongly 
because the electronegativity of the NH: group is 
lower than that of chlorine. 

Aniline has been previously studied by several 
different investigators in the ultraviolet,‘ in the 
infra-red,’ and by means of the Raman effect.® 

1H. Sponer, G. Nordheim, A. L. Sklar, and E. Teller, J. 
Chem. Phys. 7, 207 (1939). 

2A. L. Sklar, J. Chem. Phys. 7, 984 (1939). 

1941} Sponer and S. H. Wollman, J. Chem. Phys. 9, 816 
; . (a) K. Masaki, Bull. Chem. Soc. Japan 11, 346 (1936). 
(b) S. Kato and F. Someno, Sci. Pap. Inst. Phys. Chem. 
Research Tokyo, 33, 209 (1937). (c) R. Titeica, Bull. Soc. 
Roumaine Phys. 43, 15 (1942). (Abstract only available. 
See Chem. Abs. 48, 4512 (1944).) 


5 V. Williams, R. Hofstadter and R. C. Herman, J. Chem. 
Phys, 7, 802 (1939). 
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EXPERIMENTAL 


The spectrum of aniline, which lies in the 
region of 2500-3000A, (Fig. 1) was photographed 
in both the first and second orders of a spectro- 
graph of the Eagle type, using a 3-m grating ruled 
15,000 lines per inch. The light source was a 
5000-volt, 2.5-kva hydrogen discharge tube con- 
structed by the Hanovia Chemical Company. 
The plates used were Eastman Kodak 103-0 and 
103-0 uv-sensitive, the exposure times varying 
between ten minutes and two hours. The com- 
parison spectrum of iron was obtained from a 
standard arc of the Pfund type. The absorption 
cell was all quartz, 70 cm long, connected to a 
reservoir containing the aniline. After evacuation 
of the cell, the aniline pressure was controlled by 
fixing the temperature of the reservoir. Plates 
were taken with temperatures ranging between 
— 20°C and 25°C, the cell in every case being at 
room temperature. 

After each run it was necessary to carefully 
clean the absorption cell due to photodecomposi- 
tion which deposited a film on the windows. The 
film seemed to have a continuous absorption and 
so did not interfere with the development of the 
spectrum, but did necessitate longer exposures if 
allowed to build up to any appreciable thickness. 
No detectable differences were found between 
short first-order exposures and the longer second- 
order ones indicating that the photodecomposi- 
tion products did not seriously interfere. 

The aniline used in this investigation was the 


_Eastman Kodak “white label’’ grade which was 


6 (a) See reference 4(a). (b) L. Kahovec and A. W. Reitz, 
Monatsh. 69, 363 (1936). 
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Fic. 1. The absorption spectrum of aniline vapor. (Negative print.) 


redistilled in a meter-long, packed, adiabatic 
column just before use, the middle third of the 
distillate being kept for this investigation. Its 
boiling point was found to be 173.1°C at 750 mm. 

The plates were measured on a Gaertner 
.comparator, and on reduction the results from 
both orders were found to be in good agreement. 
The wave-lengths as given in Table I are the 
average of three second-order plates. Almost all 
of the stronger lines appeared as close, well- 
resolved doublets (Fig. 2) and because these were 
very sharp with no apparent degradation, the 
center of gravity of each of the components was 
measured. This method was also applied to those 
bands which did not appear as doublets. The 
precision varied from 0.3 cm for the sharp 
doublets to 2 or 3 cm for the weaker, more 
diffuse bands. 


DISCUSSION 


Over 500 bands were observed in a frequency 
range of just over 6000 cm~. In the central 
region of the absorption spectrum the average 
separation of the bands was about 8 cm~. Since 


0 +492 +292 
0+797 Geass) 


with this separation, it was impossible to make an 
unique assignment, and since the doublet spec- 
trum mentioned above stood out on the basis of 
intensity and character, an analysis was made of 
this spectrum only.’ 

From the observed frequencies of the lines, a 
difference table was prepared and the appearance 
number of the prominent differences noted. 
Several of the strongest bands were separated 
from the 0—0 band by frequencies which had 
high appearance numbers in the difference table. 
These frequencies were then taken as the vibra- 
tional frequencies for the excited state, and the 
assignment of the remaining doublets made on 
that basis (Tables I and I]). 

The first strong doublet and the most intense 
of all is the one at 34,032—34,035. Following 
Sponer and Wollman, the lower frequency com- 
ponent was chosen as the 0—0 line despite the 
fact that there was no difference in intensity 

7 Sponer and Wollman have observed similar doublets in 
chlorobenzene. Benzene also shows this doublet structure. 


The doublet separation of 3 cm™ is of the right order of 
magnitude to be attributed to rotational structure. 


0 +492 +292—42 


01948 we 


Fic. 2. The fine structure of aniline absorption bands in the region of 34,800 cm™. 
(Negative print X 4.) 
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TABLE I. The spectrum of aniline. The intensities are designated as yw= very weak, w= weak, m=medium, ms = medium 
strong, s=strong, vs=very strong, and ws =extremely strong, B=broad and D=diffuse. 








Fre- 
quency Inten- 


(em=!) sity Assignment 


Fre- 
uency Inten- 
em) sity 


Assignment 


Fre- 
quency Inten- 


(em-) sity Assignment 


Fre- 
quency Inten- 


(em~) sity Assignment 





A. The Double Bands 


33,376 w 0—231—423 


33,381 


33,501 mB 
33,506 mB 


33,569 m 
33,573 m 


33,609.1 
33,612 


33,754 
33,757 


33,801 
33,804 


33,908 
33,912 


33,948 
33,950.9 


33,980 
33,982 


33,991 8 O0—42 (953—994) 
33,994 


0—531 
0—423—42 


0—423 


34,032 0-0 
34,035 


34,094 0+62 (292—231) 
34,097 


34,193 
34,197 


34,205 
34,208 


34,241 0+292—242 
34,244 


34,283 0+292--42 
34,286 


34,297 0+797—531 
34,299 


34,314 0+282 

34,317 
0+292 
0+338 
0+797—423 
0+492—2 X42 


0+492—42 


0+492 
0+953—423 
0+492+62 
0+567 


0+717—42 





A. The Double Bands 


34,733 s 044924292 
34,736 3 —2 


34,749.2 
34,753 


34,774 ys =0-+4+-492-+4 292—42 
34,777 


34,788 0+797—42 
34,792 


34,805 0+-492+-282 
34,809 


34,814 0+-492-+292 
34,817 


34,829 0+797 
34,833 


34,861 0+-492+-338 
34,866 


34,891.1 0+797+62, 
34,895 0+292+-567 


34,923 0+-492+338-+-62 
34,926 

34,943 0+953 —42 
34,945 


34,975 0+2x492—42 
34,978 


34,985 0+953 
34,988 
34,992 


35,016 0+2x492 
35,019 


35,049 0+953+62 
35,052 


35,072 0+-797+-282—42 
35,075 


35,079 0+2X492+62, 
35,081 0+797+-292—42 


35,092 0+567+492 
35,094 


35,114 rs = 0+-797+4-282 
35,117 


35,135 0+1184—2X42 
35,138 


35,173 0+1184—42 
35,177 


0+717 


35,193 0+-953+-292 
—2X42 


35,197 


35,201 0+2X797—423 
35,205 


35,216 0+1184 
35,219 


35,223 0+953+282—42 
35,226 


35,236 0+953-+292—42 
35,239 


35,267 0+953+282 
35,270 


35,277 0+953+-292, 
35,281 0+1184+-62, 
0+797+492—42 





A. The Double Bands 


35,298 ms 0+2492+282 
35,301 ms 


35,304 ms 


35,319 0+797+-492 
35,322 


35,381 0+797+-492+62 


35,397 0+ 1446-242 
35,401 


35,437 0+1446—42 
35,441 


35,464 0+2X717 
35,467 


35,478 0+ 1446 
35,480 


35,542 0+ 1446+62 
35,545 


35,552 0+1184+-338, 
35,555 0+953-+-567 


35,626 0+2X797 
35,629 


35,690 0+2X797+62 
35,692 


35,703 0+1713—42 
35,705.4 


35,727 0+1446+4292—42 
35,730 


35,745 0+1713 
35,748 


0+2X492+-292 


35,755 0+2X797+2 X62 
35,758 


35,783 0+953+-797, 
35,786 0+1184+-567 


35,845 0+953+-797+ 62, 
35,849 0+ 1184+-567+62 


35,937 0+2X953 
35,939 j 


35,953 0+1713+4292 
35,957 —2x42 


35,993 0+1713+4+292—42 
35,996 


36,006.2 
36,009.3 


~ 04+1446+4953—423 


36,024.6 
36,027.7 


0+953+-797 
+282—42 


36,067 0+-953+-797+282 
36,069 


36,117.7 
36,120.7 


0+2X797+292 


36,221 0+ 1446+492 
+2 


36,224 92—42 


36,229 0+ 1446+797—42, 


36,231 0+2X953+292, 


0+1184+953+62, 
0+797+953 
+492—42 





A. The Double Bands 


36,260 ms 0+1446+-492+-292 
36,263.0 ms 


36,273.5 ms 
36,276 ms 


36,336 m 
36,340 m 


36,697 
36,700 


36,824 
36,828 


36,933 0+1713+ 1184 
36,936 


37,515 
37,519 


0+ 1446+797 


0+ 1446+-797+-292 


0+-1713+-953 


0+1713+-797+-282 





. The Single Bands 


33,461 wB 
33,530 m 
wB 

wBD 

mD 

w 

w 

wD 

w 

wD 

w 

w 

wB 

wD 

m 

msB 

mB 

mB 

w 


33 
3 bis 
BesesesssaSosieSes 


8es3 
we ws & 


3 
essis& 


34,448 wD 
34,458 msBD 
34,474.6 mD 
34,500 mD 
34,550 mBD 
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TABLE I.—Continued. 














Fre- Fre- Fre- Fre- 
quency Inten- : quency Inten- quency Inten- quency Inten- 
(em=!) sity Assignment (em-!) _ sity Assignment (cm-!) _ sity Assignment (em™) _ sity Assignment 
B. The Single Bands B. The Single Bands B. The Single Bands B. The Single Bands 
34,557 m 35,585  sBD 0+2X797—42 36,430. 
34,570 m 35,603 msB ” 36488 ‘ . or 270 r 
34,580  mD 35,617 ms 36,446 w 37.291 nw 
34,582 vw 35,652 ow 36,464 mB 37,304.4 mBD 
yr = pot - BD 36,477 mB 37,321 vw 
sett i ‘a a ms oe os 
oreo » = mo} ms yo med 37,369 vw 
’ - , w 36 5 ms. 37, 381 vw 
34,673  msD 35,716 i : 
34,676.8 msB 35721 sw = << 37,3808 mBD 
34,682 msB 35,741 wB 0+953+797—42, | 36568 mB 37°420.6 rw 
34,686.8 msB 0+1184—531 0-+-11844567—42 aL D 
36,581 veB 37,441.6 wBD 
34,690.7 m 35,769.6 wB 0+1446+292 36,608  mD 37 459 a 
34608 3708 — = 7468 tw 
34.719.1 mB 35,797 mm 360 On pf oe 4 
‘ ‘635. 7,532.0 msBD 
3757.6 mD as 36,045 aD 37.546 w 
crt y ' m 7, m 
Ht ae —_— = 36,672 mD 7611 
34,796 wB 0+1184—423 35,828 w 6680 37622 wBD 
34,823 msD 35,833 8B } - 7, 
34,844 wB 35 840 m one *D on w 
34,847 w 35,854 m , = , ow 
4 " 36,705.8 mD 37,671.8 vw 
34,852 ° 35,876 = maD 36.721.0 msD 37.683.6 vw 
34,857 m 35,889.4 mD ; 
34,873 ms 35,893 ms 36,735 8 37,718 vu 
34,878.6  mD 35,896.2 ms ee ome a7 6 
34,886 w 35,905.5 36,766.7 mBD 37,784.9 vu 
34,898.7 =m 359163  m 6,786 ms 37,798.3 ww 
34,002 w 35,022 ww — —_— = 
x ¥ ‘ w 37,837 vw 
34,918 vw 35,930.2 mD 36846 msBD ems we 
34,936 mBD 35,948 pe aes 4 yon ” 
35000 pl esd a 36,881 | —< 37,912.3 ww 
35,005 rw 35,981 rw 888.5 ms 37,945 vw 
35,008 vw 35,987 ww 36,905.5 m 37,952 w 
35,030 ow 36,038 mB 0+1713+292 ll — - 
35.055 m — 6947 8 w 38,023 = w 
35,126 mBD 36,089 8 36,958.8 =m 0430 bw 
35,155 m 36,092.4 8 36,976 w 38,064 w 
35,166 ms 36,096 3 oo m 38.071 vw 
35,210 1 J m 7 vu 
35931 = w e138 Os 37,018 tw 38126 ou 
35,246 w 36,157 mBD 37,081 m 38,158 u 
35,253 mBD 36,170 3B 37,048 més 38,209 vw 
35,314.1 m 36,180 ms 7,056 ow 38,221 vw 
rel we 3 2 an 37082 cma Cl 
35,341 wsB 0+9534+292+6 6, ms / m , w 
ae ue wie’ is 
A ms » 8 ’ m y vw 
35,369 w 36,245 msB 37,113 m 38,429 vw 
35,410 w 36,253 msB 37,121 mBd 38,560 w 
As m on w a mB =e vw 
35,425.0 ™m 6,294 8 . m 38, vw 
35,452 mD 36,301 w 37,151 vw 38,846 w 
35,458 m 36,310. m 37,170.9 m 38,896 vw 
35,471 w 36,318 w 7,180 m 38,944 w 
35,483 m 36,325 ms 37,187 vw 39,046 vw 
35,509 m 6,353 ms 37,196 w 39,216 vw 
35,515 ww 36,371 ms 37,202 ow 39,278 vw 
35,526 msB 0+2X717+62 36,389 wB 37,209 msD 39,341 vw 
35,563 w 36,402.1 m 37.2272 m 39,441 w 
35,575 wB 36,415 8 37,244.0 ms 39,483 w 

















between the two. Masaki has chosen 34,034 cm7 
as the 0—0 band while Titeica has assigned the 
frequency 34,035 to this transition. A strong 
diffuse band 9 cm~' to the red (34,023) appeared 
in aniline as was the case for monochlorobenzene. 
Similar satellites are associated with a number of 
the other strong doublets. 

The first doublet to the red is 42 cm from the 
0—0 frequency. This might be interpreted as a 
1—1 transition between vibrational frequencies 


994 in the ground state and 953 in the excited 


state. This 994 frequency corresponds to 991 
reported by Masaki in his ultraviolet investiga- 
tion, and to the 993 found in the Raman in- 
vestigations. Since the investigations herein re- 
ported were carried out with relatively low 
aniline pressures, the spectrum did not extend 
below 656 cm! to the red of 0—0, and hence 994 
was not found. But with 0—2X42 and 0—3 X42 
occurring in the spectrum, the non-appearance of 
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TABLE II. The vibration frequencies of aniline. (The 
letters in parenthesis indicate possible associated ground 
and excited frequencies.) 








Ground State 








Ultraviolet Symmetric Raman Excited State 
This Kahovec 
re- and This 
Masaki Titeica search Masaki Reitz research  Titeica 
_ 233 231 _ 234 62 147 
_ _ _ _ 386(a) 282 224 
421 423 423 _ _ 292 291 
529 521 531 530 534(b) 338(a) 
_ _ _ _ 615(c) 492(b) 491 
821 821 ‘ 815 817(d) 567(c) 505 
991 993 994(e) 717 
1032 1028 1030 797(d) 
1265(f) 953(e) 952 
1601(g) 1184(f) 
1446(g) 
1713 


2583 








0—994 in considerable strength weakens the 
above analysis. 

The frequencies 231 and 531 were observed by 
other investigators in the ultraviolet and also in 
the Raman effect. The frequency 423, while also 
listed by both Masaki and Titeica, does not 
appear in Raman data. 

The first doublet towards the violet is 62 cm~! 
from 0—0. This is interpreted as being a 1—1’8 
transition between 231 in the ground state and 
292 in the excited state. The frequencies 338, 492, 
567, 797, 953, 1184, and 1446 of the excited state 
may be tentatively associated with the sym- 

8 The 1—1’ transition is to be interpreted as a transition 


between one mode of vibration in the ground state to a 
different mode of vibration in the excited state. 


metric Raman lines 386, 534, 615, 817, 994, 1275, 
and 1601, respectively. The association in the 
case of the last two frequencies is less certain than 
the others since this involves rather large drops in 
frequency of 91 and 155 cm™ in going from the 
ground to the excited state. 

There remains then the frequencies 282, 292, 
717, and 1713 unexplained, 292 being also re- 
ported by Titeica. The two lower frequencies 
might be difference frequencies between unas- 
signed levels, however, these are two of the 
strongest bands in the spectrum. Similarly 717 
and 1713 might be possible combination fre- 
quencies. 

This leaves 1030 as the only symmetric Raman 
frequency that has not been associated with any 
frequency in the excited state. Since Masaki also 
found this as a 1—0 transition in the ultraviolet 
spectrum, its existence cannot be doubted. A 
possible association is the 984 frequency, which 
has here been identified as 2492. 

Several of the stronger bands which appeared 
single have been given as assignment based on 
frequencies used in the analysis of the doublet 
spectrum. 

Acknowledgment is hereby made for the special 
Graduate Research Fund grant given for con- 
struction of the spectrograph, and for supple- 
mentary grants made by the Research Institute 
of The University of Texas in support of this 
work. 
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The Thermodynamics of High-Polymer Solutions: I. The Free Energy of Mixing of 
Solvents and Polymers of Heterogeneous Distribution* 
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The theories of Flory and Huggins for the free energy of 
mixing of a homogeneous chain polymer of uniform mo- 
lecular weight with a single uniform solvent have been 
' extended to the case of a polymer mixture of varying chain 
lengths with a mixture of solvents. By making the similar 
assumptions as those of Huggins, and utilizing familiar 
statistical mechanical methods, the partial molal free 
energy of mixing of the solvent is found to be 


AFo=RT[In got (1 —Go)(1—1/mn) + u(1 — G0)? ], 


where ¢o is the volume fraction of solvent, my a simple 
function of the number average molecular weight, and uv a 
constant characteristic of the polymer-solvent mixture 
(consisting largely of a heat term, but also including y, the 
coordination number of the rubber segments). By assuming 
that a mixture of two solvents behaves like a new homo- 
geneous liquid a method of calculating u for such mixtures 
is developed. Applications of these formulas to solubility 
and fractionation are shown in a subsequent article. 





INTRODUCTION 


HE free energy relations of high-polymer 

solutions were first calculated independ- 
ently by Flory! and Huggins? for the case of a 
homogeneous chain polymer of uniform molecular 
weight in a single uniform solvent. In practice, 
however, one is never dealing with a polymer 
sample in which all the molecules have the same 
molecular weight, but rather with a mixture of 
molecular weights, forming a virtually continu- 
ous distribution. If, on the other hand, one 
intends to extend the thermodynamic treatment 
of such solutions to the problems of fractionation 
and precipitation it is of interest to extend the 
treatment to a solvent medium consisting of a 
mixture of two pure solvents. The purpose of this 
and the following paper is to extend the theory to 
these two cases and to give a thermodynamic 
treatment of solubility, fractionation, and pre- 
cipitation. Although the final results of Flory and 
Huggins are fairly analogous, we preferred to 
base our derivations upon the development of 
Huggins who, in our opinion, makes less, and 
more justifiable, assumptions concerning the 
nature of the rubber molecule than does Flory. 


1. Entropy of Mixing of Chains of Varying Length 
with a Solvent 


Suppose a mixture consisting of ; molecules 
of molecular weight m (in units of arbitrary 


* The major portion of the following work was done on a 
grant from the Rubber Reserve Company for fundamental 
research in the field of synthetic rubber. 

1P. J. Flory, J. Chem. Phys. 10, 51 (1942). 

2M. L. Huggins, Ann. N. Y. Acad. Sci. 43, 1 (1942). 


submolecules, the density of all molecular weight 
fractions being taken as the same),* mz molecules 
of molecular weight m2, 13 of m3, etc. Then the 
total number of sites, for submolecules and 
solvent molecules is 


N=Notnim+nome+---+nymy, (1) 


where N)=the number of sites occupied by the 
solvent molecules. 

In order to calculate the entropy of this system 
we have to determine the total number of 
configurations. We shall proceed in the same way 
as Huggins did for the binary system.‘ 

Suppose now that ; molecules of the type m, 
n2 of the type mz, etc., up to 2;_1 of the type m;_; 
have already been put into the “‘lattice’’ as well 
as K—1 molecules of the type m;. 

Let vix,, be the number of alternative sites 
available for the /Jth segment of the Kth molecule 
of the type m;. Then for the first segment of the 
Kth molecule of the m; type 


vi,x,1= N(1—f)), (2) 


5’ Each submolecule is considered as the same size as a 
molecule of solvent so that m; is in reality the ratio of the 
molecular volume of the polymer molecule of molecular 
weight M; and the molecular volume of the solvent Vo or 
mi = M;/pVo. lf p, the density of the polymer, is considered 
as independent of the chain length, then m; is, for any given 
solvent, diréctly proportional to and therefore a measure 
of the molecular weight. 

4 There is, however, one difference between our treatment 
and that of Huggins. Huggins considers each molecule of 
molecular volume m as consisting of / freely orienting 
segments each of them composed of p submolecules. Since, 
however, in the final formulation / and p appear only as the 
product /p=m, we did not use this decomposition and con- 
sider simply the molecule as composed of m segments. 
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where f;=the fraction of the total number of 
sites already occupied. 


i—1 
> nym;+(K —1)m; 
i=1 


/= ; 3 
f = (3) 





Once the first submolecule is placed in a site So, 
the position of the other submolecule is greatly 
limited. Huggins points out that if the immedi- 
ately adjacent site S; is occupied, two of the sites 
immediately adjacent to S; must also be occupied, 
provided S; is not the end of the chain.5 Hence 
the probability that both S; and So are occupied 
by adjacent molecules of a previous chain is: 


P=2f;/7, (4) 


where y is the ‘coordination number,” i.e., the 





number of neighbor sites around each site. But in 
our case Sp cannot have been occupied by a 
submolecule of a previous chain since we assumed 
that it is occupied by a first submolecule of the 
present chain. Hence the probability that 5S, is 
unoccupied and free for the second submolecule 
is greater than (4) and is 


P=(1-f;)/(1—(2f;/)], (5) 


and the number of sites available for the second 
submolecule is 


vi,x,2=7(1—f))/[1—(2fi/y)], (6) 
and similarly for the third submolecule 
1n.x,3=(¥-1)(1-f)/T1-Q2fi/)], 


and so on. 


Hence for the Kth chain of the ith type the total number of distinguishable configurations is 


vi,K=(N/o)(1—fi)v(y—1)"" [1 —f5J/[1 — (2fi/) I, (8) 


where o is the symmetry number of the chain (i.e., 


the definition of f; (3) 


(N- 


vi.k=¥(y—1)™*- 


2 if the ends are interchangeable), or introducing 


1 
> nym;—(K—1)m,)" 
1 


= 





o 


(N—2/yLX mams—(K 1) 


Then, the total number of configurations for the 1; molecules of the ith type may be written as: 


ni-l én% 


I] (V- nm;—km,)” 
i=1 








ni aaa k=0 
ni-l 


I] .=*x 
k=1 o 


k=0 


(9) 


Tl (N= 2/70 5 namitkm,J]yn 


i=1 


Taking the logarithm and replacing the summations by integrals (equivalent to the approximation of 


Sterling’s formula), 


“ ei ies im} ; 
In T[ »% x=; In { ————— +m. f In (N—> nym;—Km,)dK 
o 0 i=1 


K=1 


= (m1) f In (N—2/7E nam,— Km) aK. (10) 


i=1 


5 For long molecules, as Huggins has shown, this possibility introduces only a negligible correction. 








aa 
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Integrating 
4 Ai cai i i 
In [] x, xn=n; ln ( ————— } —nm;-(N-Y nm, In (N—-Y¥ nm,) 
K=1 o i=1 i=1 
i—1 (m i-1 


+(N-—D nm; ) In (V— a nm;)+ yL(N —2/y 3 nm) In (N—2/y ¥ nmi) 


i=1 mM; i=1 i=1 





i—1 i—1 
—(N—2/y > nm) In (N—2/y ¥ nym;)+2/ynm; |. (11) 
i=1 i=0 
The total number of distinguishable configurations for the whole system is 


W= CTE Tc) CT (me) JCI TD vex CTT (12) 


and the entropy of the system 
S/k=|In W=>0 In [Tx «-—D ny In 2;4+Dd ny. (13) 
a K=1 i 


Substituting (11) in (13) and simplifying 


a(y—1)"** 


o 


S/k= E min|™ |- LX niin ni/N—No ln No/N+y7/2(N—2/y L nimi) 





L nim; 
xin( 1-27 p= )-1a5cwm JL(1—2/yN- non) In (1—2/yN- _— 


i=1 i=1 


—(1-—2/yN:- a nm;) In (1—2/yN- > nm,;) |. (14)® 
i=1 

So far no new assumption or approximation has been made beyond the ones used in Huggins’”’ 

theory. However, in order to go further than Eq. (14) one has to make the following simplification : 

the last term of the Eq. (14) is very small since it involves essentially a difference of sums extended 

once to the ith molecule type, once to the (¢—1)th type. If we now assume that any given molecular 

weight fraction is small, mm (()/ nm; it is legitimate to expand the logarithms in series and cancel 
all but the first term. Then this term becomes 


—y/2:- X N/m; [(1—2/yN- mm) In (1—2/yN- ¥ mm) 


—(1-—2/yN- z nm;) In (1—2/yN- E nmS+¥ n;. (15) 


Equation (14) has taken then the simplified form: 


1 ae | alg 
S/k2Y n; In | —————— 1-0 In n;/ N— No In No/N 
i Co ry 
+y7/2(N —2/y p nm;) In (N—2/y ) ¥ nimi) +> n; (16) 
6 If one now assumes that only molecules of one chain length are present, 1; = N,, n2=n3;= --- =0 then Eq. (14) reduces 
to: 





v(y—1)™ N, No M ‘ii: ~~ No+(1—2/y)mN, 
S/k=NrI p pow’) — Ne In Eo — No ln a + (m—1)E4/2m- Not (1/2 1)N-]In | wee | 


which is identical with the equation derived by Huggins. 
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The entropy of mixing is of course 
ASmix = S—So—D Sis (17) 


where Sp and S; are the entropies of pure solvent and pure rubber fraction respectively. Substituting 
in Eq. (16) the conditions 7140 and m2=n3=m4= --- =0; 240 and m;=n3=m=---=0, etc., and 
No=0 separately one finds 


So=0, (18) 
S;/k=n; In [y(y—1)"*? )]/o+n, In m+ (y/2—1)nm; In (1-2/7) +03. (19)? 
After subtraction one obtains the entropy of mixing. 


ASnix/k = we * nN; In nim;/N—No In No/N 
++7/2(N—2/y X nim) In (N—2/y Y nimi) —(y/2—1) In (1—2/y)- 20 mim; (20) 


After partial differentiation and some simplifications, Eq. (20) yields the partial molal entropies. 














ASo 0(ASmix/R) nN; No No Y 2 2 
— = ————=4+ 2 —-In 4 (1-— 5 ns) +— 2 nm) (21a) 
R ONo iN N N 2 yN i yN i 
AS; 0(ASmix/R) n; nm; m:No 
— == +m; 5 ——In ———-1+ 
R On; iN N N 
( 2 
{1 ——— Te nym; > nymN; 
Y yN i i 
+(2-1) n ~m\ 1- ) (21b) 
2 2 N 
Pao 
\ Y y, 








If we introduce now measurable values, e.g., volume fractions, one has the relation 


No/N=¢0, nimi/N=¢i, YL nmi/N=L ¢i=1— 90. 


Further, the number average molecular weight (in units of submolecules) is: 
mn = nmi/Y mi=(L $:)/LX (bi/m:) J (22) 


and hence 
) > n;/N= 1 — do/Mwn. (23) 


Introducing this definition in Eq. (21), expanding the logarithms in the last term and neglecting 
the terms in ¢* or higher, one obtains finally: 


ASo= —R[In go+(1 —1/mwy) (1 —0) +1/7(1—90)?]; (24) 
AS;= — R[In $:+1—migpo— (m;i/my) (1 —$0) + (m:i/7) 07]. (25) 





Several remarks can be made concerning Eqs. as would follow from the calculation given above. 
(24) and (25). (1) It is not apparent why the last This substitution is made by analogy to Huggins, 
term in Eq. (25) is m;/y-¢0? and not mi¢o?/(y—2) who gives no justification for it. Actually, the 
es rie denominator has to be y and not y —2 if Eqs. (24) 


7 This is what Fl lls the disorientation entropy of p . . 
a eee and (25) are to be consistent and satisfy the Gibbs 


the polymer mixture. 
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relation. Presumably this substitution corrects 
for the fact that the expansions in Eqs. (21a) and 
(21b) are not exactly equivalent, the terms in 
(21a) being all positive while the terms in (21b) 
are alternately positive and negative. Hence when 
the higher terms are neglected, the same degree 
of approximation is not reached. (2) For a binary 
system m;=m,=m and this reduces to the 
familiar equations of Huggins: 


ASo= —R[In go+(1—1/m)o,+(1/y)o7]; (26a) 
AS,= —R[In ¢,—(m—1)¢o+(m/y)o?].  (26b) 


(3) The partial molal entropy of the solvent does 
not depend on the molecular weight distribution 
of the high polymer but only on the number 
average molecular weight my, while the partial 
molal entropy of the polymer depends on both 
my and the molecular weight distribution. This 
is important because it shows that the colligative 
properties (osmotic pressure, vapor pressure 
lowering, freezing point lowering) are independent 
of the actual molecular weight distribution of the 
polymer and always permit determination of the 
number average molecular weight. 


2. Heat of Mixing 


No exact formula exists for the time being for 
the computation of the partial molal heats of 
mixing. However, at least in the case of non-polar 
polymers, it was shown by the experiments of 
Gee and Treloar§ and by some unpublished 
experiments of the authors that the total heat of 
‘mixing can be fairly well represented by a 
Scatchard-Hildebrand formula, i.e., 


AAnix= Vin(ao— ar)*hodr, (27a) 


where a and a, are the square roots of the 
cohesive energies (internal heats of vaporization) 
per unit volume : 


ag= (Eoceony/ V,)t= (Ao—RT/ Vo)}, (28) 
ar= (Econ) / V,)i= (A,-—RT/ V,)}, 


and V,, the molecular volume of the mixture 


Vn=Xo VotX, V,= Vo V,/ (0 V.+¢- Vo). (29) 


8 G. Gee and L. R. G. Treloar, Trans. Faraday Soc. 38, 
147 (1942), 
® The subscript 0 represents the solvent, and 7 the rubber 


(polymer). 


In order to expand the Scatchard-Hildebrand 
formula to mixtures of polymers and solvents, let 
us recall at this point the basic assumptions from 
which the formula (27) is derived. The cohesive 
energy of a mixture is supposed to be the sum of 
three types of interaction, writing each in terms 
of the volume fraction 


(Econ) mix = (A 11617 +2A 126162+A 2262”) Vin, (30) 


where Aj;, Ai2 and Ag are constants charac- 
teristic of the interaction between the molecules, 
A, and A22 between molecules of the same,type, 
Ai2 between molecules of different types. 

For the pure components, 


(Eeon)1=Ai1Vi; (Eeoh)2=A22V2; (31) 


=a; Vi; =a Vo. 


For the interaction of two different type ° 


molecules, itis now assumed 
Ay= (A 311A 22). (32) 


This assumption was first made by van der 
Waals. The quantum mechanical treatment’? of 
the intermolecular forces of non-electrostatical 
nature corrected this formula to 


A325 (A11A22)}. (33) 


Assuming Eq. (32) to be valid, Eq. (30) can 
hence be written 


(Ecob) mix = [(A11)*61+(A22)*p2) ?, (34) 


which, after subsequent transformation, under 
the additional assumption that the mixture is 
ideal in the sense that volumes of the components 
are additive, leads to Eq. (27). 

In the case when one of the components is 
polar the relation of Eq. (32) can no longer be 
expected to hold even approximately. Neverthe- 
less, experiments have shown that, in this case, 
the total heat of mixing can be roughly repre- 
sented by 


AA mix = K Vin(ao— a)’ hodr, (27b) 


where K is a constant somewhat larger than unity 
and correcting for the fact that always in this 
case Ai2<(Ai1A22)!. For mixtures of non-polar 
compounds K = 1. 

It must be emphasized that the cohesive energy 
is essentially a localized affair determined by the 


10F. London, Zeits. f. physik. Chemie B11, 222 (1930). 
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interactions between adjacent atoms and groups 
of atoms belonging to the several components. It 
is essentially the interaction between the sub- 
molecules of the polymer and the molecules of the 
solvent, and is independent of the actual chain 
length. 

Let us apply now Eq. (27a) or (27b) to the 
complex system of a mixture of a solvent and of 
polymer molecules of differing chain length. 
From the outset we can assume that a, is the 
same for all the different rubber molecules and 
that therefore there is no heat of mixing from 
this source. The only heat of mixing developed is 
that between the solvent and the polymer which 
we may treat for our purposes as one component. 

Introducing the molar fractions, Eq. (27a) 
becomes 


Anix= (Xo Voth X mM; Vo) (a9—a,)? 
Xo Vo(du X m;Vo) 





x 
(XoVo+L X m;Vo)? 
Xo > Xm; 
(Xo+X Ximi) 


= Vo(ao— a,)? 





(35) 


The partial molar heats of mixing are obtained 
by differentiation of Eq. (35). 


— OH mix 
ane 
0Xo 


= Vo(ao—a,)? 





(L Ximi)? 
(Xot+>d X m;)* 
= Vo(ao—a,)*(1—¢0)?; (36) 








| compe 0A nix 
aH.=( ) 
OX ; 


: mi(X 9)? 
(XYotl X m;)?* 
= mV o(ao— a,)*oo?. (37) 





= Vo(ao—a, 


In order to be able to treat the solubility and 
other properties of a polymer in a mixture of two 
solvents, we have to extend the formula for the 
heat of mixing in order to include this case. We 


found that solvent mixtures can be treated with a 
fair degree of accuracy assuming the mixture of 
the two solvents to behave as a homogeneous 
liquid," with respect to the polymer, this new 
“liquid” having characteristic values of a and V, 
am and Vm. The values of a» and V» can be 
obtained from the theory of heat of mixing 
developed above (Eq. 34). 


(Baada/ Vm — Gins* _ ((A 11) git (A 22) 12) ’ (38) 
or 
Om = $101 + $202 = aitg2(a2—a), (39) 
while ‘ 
Vin= Vi V2/d1Veto2Vi1 
= Vi[(1+¢2V2/(Vi—V2)) J. (40) 


3. Partial Molal Free Energy 


We are now able to express the partial molal 
free energies of the system considered : a mixture 
of polymer molecules of different molecular 
weight and of one or two solvents. 

Combining Eqs. (24) and (25) with (36) and 
(37) respectively, we obtain 


AF, =RT[In ¢0+(1—¢0)(1—1/my) 
+nu(1 — 0)” ], (41) 
AF;=RT[In ¢:+1—m ido 
—(m;/mn)(1—0) +m indo? ], (42) 
where p= 1/7y+K Vo(ao—a,)?/RT in the case of a 
simple solvent and 
KV;(a:—a,)? 


RT 





bm = 1/y+ 


ade—-Qj 2 Vi- Ve 
| 1+ )| /|1+6 | (43) 
Qi1—-a, Vo 


in the case of a mixture of solvents (1) and (2). 

In the following paper these expressions will be 
applied to the problems of solubility and fraction- 
ation of a rubber-like polymer. 

The authors wish to thank Professor Hugh S. 
Taylor for his constant encouragement and 
advice, and the Rubber Reserve Company for 
permission to publish this material. 








1 This can be true only if the interaction of the solvents 
with the polymer is of second order in comparison with the 
interaction between the two solvents. 
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The free energy relations for heterogeneous molecular weight distributions developed in 
Part I, are applied to problems of solubility and fractionation. Critical conditions for solubility 
are obtained. A rigorous expression for the solubility is derived, although certain approximations 
are made to facilitate calculation, and to permit extension to polymers which are part “gel.” 
Following the same methods the thermodynamic equilibria involved in fractionation are 
described by mathematical expressions which permit a comparison of the extraction and pre- 
cipitation methods. The effectiveness of fractionation is shown to be strongly dependent upon 


concentration. 





INTRODUCTION 


LORY! and Huggins? first developed theories 
of the entropy of mixing of polymers of 
uniform molecular weight with solvents and soon 
thereafter* 4 extended their treatment to include 
the solubility of such homogeneous polymers. 
Later, using a modified form of the Flory theory, 
Gee® attempted to give a rough treatment of 
fractionation. Until now, however, no systematic 
treatment has been given for polymers consisting 
of a heterogeneous distribution of chain lengths. 
Such distributions are to be found in both natural 
and synthetic rubbers and in all similar chain 
polymers. 

It has been shown in Part I* that in a mixture 
of solvent with such a heterogeneous polymer, 
the partial molal free energies of mixing of the 
solvent AF» and of the ith fraction of the polymer 


AF; may be represented as: 


1 
BFy= RI In bo (1-4) (1-— 


7 
+n(1 -+4)'| (1) 


* The major portion of the following work was done on a 
grant from the Rubber Reserve Company for fundamental 
research in the field of synthetic rubber. . 

1P, J. Flory, J. Chem. Phys. 9, 660 (1941). 

2M. L. Huggins, J. Chem. Phys. 9, 440 (1941). 

3P. J. Flory, J. Chem. Phys. 10, 51 (1942). 

4M. L. Huggins, J. Phys. Chem. 46, 151 (1942); Ann. 
N. Y. Acad. Sci. 43, 1 (1942). 

5G. Gee, Trans. Faraday Soc. 38, 276 (1942). 

6R. L. Scott and M. Magat, J. Chem. Phys. 13, 172 
(1945). 
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mM; 
ofl — $0) 
N 


m 


BF.=R1|In 6.+1— 
— mibo+ umd’ | (2) 


where ¢o and ¢; are the volume fractions of 
solvent and ith fraction respectively. 

m; is the ratio of the molecular volumes of the 
ith fraction V; and that of the solvent Vo and, 
assuming that the density of the polymers is 
independent of the chain length, a direct measure 
of the molecular weight. 


m;= V;i/Vo= Mi/p, Vo. (3) 
Mn is the number average m, defined as 
Le nim; } 2 mM; 
— : 0 





on Lmi/d: 


and related to the number average molecular 
weight by a relation similar to Eq. (3). 

uw is a constant for any given system, consisting 
of the heat of mixing and a term from the entropy, 
such that 


w=1/y+KVo(a0—a,)?/RT, (5) 


where y=the coordination number of the seg- 
ments (about 3-4), ao, a-=constants for the 
rubber (polymer) and the solvent (the square 
root of the internal energy of vaporization per cc) 
and K an empirical constant relating to the 
polarity of the substances in the mixture.” 


7 For a more comprehensive discussion of these constants, 
see Part I. 
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Equation (1) is identical with that for the free 
energy of the solvent in a binary system as 
derived by Huggins, with the exception that my 
replaces m. Equation (2), however, is not identical 
with that of the rubber in the homogeneous 
system. 


A. Solubility of a Polymer of Distribution f(m): 
Distribution of a Fraction of Molecular 
Weight m; Between Solution 
and Swollen Precipitate 


From Eqs. (1) and (2), one may deduce the 
solubility conditions and the distribution of the 
rubber between two phases. 

If one plots AF against decreasing ¢o, one finds 
that for sufficiently small values of yn, AF, de- 
creases uniformly, signifying complete miscibility 
of solvent and polymer. (Fig. 1.)§ With a larger 
value of u (that is, a larger heat of mixing), how- 
ever, the curve shows a minimum and a maximum 
indicating only partial miscibility. There will 
therefore be a critical value of yu, u., at which the 
solution breaks into two phases. For this value 
the curve will show a point of inflection, and 
applying the familiar conditions that, here, the 
first and second derivatives of AF» with respect to 
¢o must vanish, one finds: 


e= (1/7) +(K Vo(Aa)*/RT) 
=1/2[1+(1/my)*}?. (6) 


This is identical with the similar equation,* 4 
for a homogeneous polymer except that my 
replaces m. In short, all the critical conditions are 
the same as those for a binary system if one 
utilizes the number average molecular weight. 
These critical conditions depend only on my and 
not on the precise nature of the distribution 
function, which may correspond either to a fairly 
narrow, homogeneous fraction, or to a very broad 
distribution. 

The critical point corresponds also to a critical 
value of ¢» and of ¢, (the volume fraction of the 
total rubber.) 


8For identical curves for polymers of homogeneous 
distribution, see M. L. Huggins, Cellulose and Cellulose 
Derivatives, edited by Ott (Interscience Publishers, Inc., 
New York), p. 893. 





mn 
We= = ; 7 
™ 1+ /my 1+(1/+/my) - 
rje= Ce. 7b 
($r) (2 $4) re (7b) 


Similarly one may speak of a critical solution 
temperature 7, or a critical value of (ao—a,), 
(Aa)c. 


T.= {K Vo(Aa)?}/{R[1/2(1+(1/+/mw))? 
—(1/y) I} (8) 
(Aa).= { RT[1/2(1+(1/+/my))? 
—(1/y)}}/KVo. (9) 


Certain consequences follow immediately from 
these equations: (1) As pointed out above, the. 
critical point is dependent only upon the number 
average molecular weight. Since, for values of u 
less than yu., no separation of a rubber phase is 
possible, it is clear that the presence of low 
molecular weight material in the polymer tends 
to keep the high molecular weight material in 
solution by decreasing the value of my. In short, 
low molecular weight material actually makes the 
higher fractions more soluble. This factor is crucial 





+0.4 . 7 * pee T 


+0.2} 7 








AF, (calories/mole) 














L 1 4 l 
0 0.10 0.20 0.30 
or 
Fic. 1. Partial molal free energy of solvent AF as a 


function of uw and @¢,;, the volume fraction of rubber. 
mn = 100 (my=10,000). 
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in any consideration of fractionation or precipi- 
tation curves. (2) Since yu, is inversely pro- 
portional to T (Eq. (5)), the solubility increases 
with an increase in temperature. (3) If for a given 
rubber-solvent system, (ao—a,) is less than the 
critical (Aa), in Eq. (9), the rubber and solvent 
are necessarily completely miscible in all pro- 
portions. 

To determine the solubility at any other point 
than at the critical point, one must utilize the 
equilibrium conditions between the two phases: 


AFv=AFy’, (10) 
AF;=AF/, (11) 
where the prime represents the swollen rubber 


phase. 

There are, of course, 7 simultaneous equations 
of the type of Eq. (11), one for each distinct 
fraction. 

In other words 








1—do 
In do+1—d0— +pu(1—¢o)? 
mn , ; 
=In go’ +140’ -——+u(1-0')*;_ (12) 
1N 


mM; 
In ¢;+1——(1—¢0) —mido+mindo® 


UL 
mM; 
=In $+ 1———(1 — $0’) —mido' +mindo”. (13) 
mn 


Note that the distribution functions in the two 
phases are presumed to be different ; hence my’ is 
not the same as my. 

Solving Eq. (13) for In ¢,’/¢; and substituting 
from Eq. (12), one obtains 


pi go 

In —=mz{ 2pn(¢0—¢0) —In i (14) 
i go 

or 


¢;'/¢:=exp (Ami), (15) 


A =2u(0— 0) — (In G0/¢0’). 
A relation similar to Eq. (14) is given by Schulz :° 
(16) 


where 


C= Ke-P (at By)! 6, 
where C=solubility of material of molecular 


*G. V. Schulz, Zeits. f. physik. Chemie A179, 321 (1937). 
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weight M; P=degree of polymerization of 
polymer = M/Mo; Mo= molecular weight of mon- 
omer ; y =concentration of precipitant (not the y 
of Part I of this paper) ; 8=a temperature func- 
tion; A, B, K=constants. 

This equation yields the same exponential de- 
pendence of the solubility on the molecular 
weight as in Eq. (15), but in our opinion, the 
assumptions made in the derivation are faulty, 
and the significance of the constants uncertain. 
Schulz’s work is important, however, as the first 
attempt to treat the problem of solubility and 
fractionation from a thermodynamic standpoint. 

Equation (15) represents the fundamental 
basis of the theory of solubility and fractionation 
presented in the following pages. For any given 
system, A is a constant; and, therefore, Eq. (15) 
may be used as a strictly valid relation, showing 
the distribution of material of different molecular 
weights between the solution and the rubber 
phases. 

We may now undertake to express the solu- 
bility. Let the actual distribution of molecular 
weights be represented by a distribution function : 


dw;= f(m)dm,, (17) 


where dw; is the fraction!® of the total rubber of 
molecular weight m; and f(m) is a normalized 

function, such that 
f f(m)dm=1. (18) 

0 

What is the solubility if a certain volume 2, of 
this rubber is placed in a volume vp of a certain 
solvent? The total amount of rubber, v; in any 


given fraction is the sum of the amounts in the 
solution phase and in the swollen rubber phase. 


Vi=PMs+Gi Up, (19) 


where v;=v,dw,;=v,f(m)dm;; v,= volume of solu- 
tion phase ; v,= volume of swollen rubber phase. 
Substituting the condition of Eq. (15), 


Vi=ois+ow, exp (Am); (20) 
bi=0;/(Vst+v,e4"1) =0,f(m)dmi/(vs+Vp exp (409s 


10 Either weight fraction or volume fraction, since the 
density of the rubber is presumed to be uniform, inde- 
pendent of molecular wae. This assumption however, is 
only approximate. It means in reality that we neglect the 
end atoms when compared to the length of the chain. For 
high molecular weights, this is obviously justified. 
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Fic. 2. The approximation constant A. as a function of yp. 





0.75 0.80 


The total solubility in terms of volume fraction 
of rubber is therefore 


Pom, en J v,f(m)dm/(v.-+0,e4"). (22) 


Equation (22), while rigorously correct is not 
very helpful when it comes to actual calculation. 
Given the distribution function f(m), either ex- 
plicitly or graphically, the ratio of v, to vp and the 
value of uw for the mixture, one can make the 
calculation, for v, and v, are functions of the 
above variables—but the computation may be 
exceedingly long and tedious. 

Exact calculation of A is extremely difficult, 
since both ¢» and ¢»’ are functions not only of my 
and my’ but also of the distribution function 
itself. It is to be noted, however, that the values 
of ¢, (i.e., 1—@o’), calculated for various m’s, 
approach the limiting curve for m= asymp- 
totically as uw increases. Hence, for values of yu 
sufficiently above the critical point, one may take 
for A, the value A. calculated for m=. For 
m=, go=1, since the rubber is completely 
insoluble above the critical point, so that 


Aw = 1($r') mae +In [1 - (br deren} (23) 





1000 T T 


100 


Bo 


10 











0.50 0.55 0.60 0.65 0.70 0.75 0.80 
7 


Fic. 3. The approximation constant 8.. as a function of u. 
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Since, when m= © in Eq. (1) (¢,’) mes is a func- 
tion of u alone: 


In [1— (6+!) mao J+(62’)mantit(r’)mmo=0, (24) 


one may calculate A. as a function of yu. The 
results of these calculations are shown graphically 
in Fig. 2. Actually, these calculations of A are 
accurate virtually to the critical point, thanks to 
compensating errors. 

Certain other approximations may also be 
made: For dilute solutions, where the swollen 
rubber phase is small in comparison to the solu- 
tion, v, may be taken as approximately equal to 
vo, the amount of solvent added. 

Further, one may write v, as a function of the 
amount of rubber in the precipitate, v,’ and the 
percentage of the precipitate which is rubber, ¢,’. 


Vp=0r /br =v,'/(1—dv'), (25) 


and hence the solubility 





or= 


c) J d 
f v,f(m)dm (26) 


vot(v,'/¢,’) exp (Am) 


We may, however, determine an approximate 
limiting value of ¢,’ by taking the value for 











My =, which is the asymptote for all finite 
values of my as uw increases. The asymptotic value 
1/(¢,’)m=0, Which we shall henceforth call 8. is a 
function of u only, and may be calculated from 
Eq. (24). ; 

Values of 8. are shown graphically in Fig. 3. 

A further approximation may be made in cases 
where yp is appreciably greater than the critical y,. 
Then the solubility will be small and nearly all 
the rubber will be in the precipitate, and one may 
write 


Nearer the critical point, the solubility is con- 
siderable, and v,’ is appreciably smaller than »,. 
Here also, however, ¢,’ is less than (¢,’)m=., and 
partially compensates. Since this deviation of ¢, 
is not usually as marked" as that of v,, one may 
say as a general rule 


Bo,’ =v, /o,’ S0,/(br’) mao = Bao r- (28) 


Therefore one may write as a final approximate 
form for the solubility 


“ f(m)dm 
a . 29 
. J (v0/v,) +Bx exp (A om) ' 


Because of the condition of Eq. (28), one will 
expect that the actual solubility will be greater 
than that calculated from Eq. (29). Errors in 
determining A are distinctly second order com- 
pared to those involving 8. 





B. Solubility of ‘Gel’? Rubbers 


One further point must be noted—Eq. (29) 
as it stands applies only to rubber which is 100 
percent soluble in ‘‘good’”’ solvents (u<0.50)— 
the so-called ‘‘sol rubber.’’ Actually much rubber 
consists of more or less highly cross-linked ma- 
terial which is completely insoluble in all solvents 
—the so-called ‘gel rubber’—all vulcanized 
rubber is of this type. This gel can swell greatly 
so that the rubber phase may be very voluminous, 
but none of the cross-linked material will be found 
in the solution phase. 

This swelling may be calculated, as was shown 
by Flory and Rehner.” By assuming the cross- 

11 [t can be shown mathematically that, even at the crit- 
ical point, ¢,’ is greater than }(¢r’) moo. Hence Bo <8 < 3B. 


122 P, J. Flory and J. Rehner, J. Chem. Phys. 11, 521 
(1943). 


V,/20,. (27) 
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linked gel to be an infinite three-dimensional 
network of chains and deriving a suitable’ ex- 
pression for the entropy of elastic stretching, they 
obtain, as the,condition for swelling equilibrium 


Pr Voo,t 





In (1 — or) +or+ugr?+ =0, (30) 


c 


where V)= molecular volume of solvent ; p,=den- 
sity of rubber; 1/,=average molecular weight of 
chains between cross links. 

For gels which are not highly cross-linked, M, 
is very large (of the order of 100,000) and the last 
term of Eq. (30) is negligibly small. If this last 
term is neglected, the equation is identical with 
Eq. (18). Hence for loosely cross-linked networks, 
such as found in unvulcanized rubbers, the degree 
of swelling may be taken as roughly the same as 
that of the sol for the range in which we are 
interested .% 

We may therefore write as the generalized 


case: Kl \d 
- m)am 
nok . 31 
i J (vo/v,) +Be4™ ~ 


where K is the fraction of the rubber which is 
*‘sol.”’ 

The solubility may also be expressed as the 
fraction of the rubber dissolved, S. In this case 
one may write 











ors 
S= fraction dissolved = 
Ur 
” m)dm 
=<=K i) A . (32) 
o 1+(v,/v0)Be4™ 


C. Fractionation 


(1) Methods of Fractionation 


For a long time, at least two different methods 
have been in general use for fractionation of high 
polymers : 

(a) Precipitation method.—The raw polymer is 
dissolved in a good solvent, and then successive 
amounts of a non-solvent are added to the solu- 
tion, producing fractional precipitation. The 
highest molecular weights are the first to precipi- 
tate, then the next highest, and so forth. 


13 For lower values of u (below 0.50), the sol is completely 
soluble, but the gel is never soluble. 
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(b) Extraction method——The raw polymer is 
successively extracted with mixtures of a good 
solvent and a non-solvent, increasing the relative 
proportion of solvent with each extraction. In 
this case, the lowest molecular weights are ex- 
tracted first, and then higher molecular weights 
and so on. 

Both of these methods can be interpreted in 
terms of the same theory, based essentially on 
Eq. (15). In both cases, each step consists of the 
establishment of thermodynamic equilibrium be- 
tween a swollen rubber phase and a solution con- 
taining varying quantities of solvent and 
non-solvent. 


(2) Theory of Fractionation 


In a polymer of a given distribution f(m) 
(Eq. (19)), the equilibrium distribution of any 
fraction of molecular weight m between solution 
and precipitate phase is expressed by Eq. (15). 

We may then develop equations in a similar 
manner as before. The actual amounts of the ith 
fraction in each phase w; (solution) and w,’ may 
be written as 


Vi=$0.+6;/0,=U;+Uy ; (33) 
wi! /wi=$;'Vp/bis=(vp/v.) exp (Am). (34) 


Solving Eqs. (33) and (34) for w; and w,’ and 
utilizing the approximations of Eq. (29) 


Vv; 
/ 


w,' = 
1+(v./v,) exp (—Am,) 





Vi 
~ 1+(vo/0,'B) exp (—Am,)’ 
U; 


WwW;= 
1+(vp/v.) exp (Ami) 





(35) 





= “ (36) 
1+(v,’/v0)B exp (Am;,) 


In all our further considerations, we shall 
assume that the solvent-non-solvent mixture may 
be considered as a homogeneous liquid, the same 
in both phases,'‘ and with an interaction with the 





4 This is not strictly true, for the solution will actually 
be slightly richer in non-solvent and the precipitate in 
solvent, but treatment of such a condition leads to tre- 
mendous complications, and the above assumptions seem 
to introduce no errors of great significance, when compared 
with unpublished experimental results of the author. 


183 


rubber characterized by a yum calculated from 
Eq. (44) of Part I. 

From this um, which depends on the relative 
proportions of solvent and non-solvent, one may 
calculate A and £ as before. v,’/vo is nothing more 
than the concentration c (in units of volume 
fraction) of the precipitated rubber in the total 
mixture. 

We may therefore express the distribution 
functions of the rubber in the solution and in the 
precipitate as follows: 


Precipitate : 


f(ms) =Cf(mi) ]/1+Lexp (— Am) }/Bc, 


Solution : 
f.(m;) =(f(m;) /(1+ 8c exp (Am;) ], 
where of course: 
fp(mi) + f.(mi) = f(mi). 


The actual amounts in precipitate and solution 
are therefore: 


(37) 


(38) 


(39) 


* — f(m)dm 
wan 1+ (e-4")/Bc 


f f(m)dm 
W,=U, a er eae 
0 1+8ce4™ 





(40) 


(41) 


It is apparent immediately'® from the above 
equation that in any such distribution of the 
rubber between the two phases, most of the high 
molecular weight material will be in the precipi- 
tate, and most of the low molecular weight 
material in the solution. There is no sharp divi- 
sion since there will be some of every fraction in 
both phases. In the case of the high molecular 
weights, the amount in the solution may well be 
negligible since 


lim 1/(1+8ce4™) =0, 


ma 


(42) 


but there will always be at least Bc of the low 
molecular weight material in the precipitate, 


16 All of the following qualitative conclusions are quite 
rigorous. Actually Eqs. (40) and (41) are precise until one 
makes simplifying assumptions in order to calculate 8 and 
A. The dependence on molecular weight is quite accurate, 
however. 
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since 


1 Bc m 
m»o 1+(e-4"/Bc) 1+ Bc v,+2, 


This last fact is important, for in the precipi- 
tation method of fractionation, a finite amount of 
the low molecular weight material will always be 
in the precipitate, a fact which is well known 
empirically. That this should be is to be seen from 
a consideration of the interaction between the 
molecules. If the interaction between the rubber 
molecule and the solvent is the same as the 
interaction between it and the rest of the rubber, 
then it is completely indifferent as to its distri- 
bution, in which case, its concentration will be the 
same in both phases. In general however, the 
rubber molecule interacts more strongly with 
molecules of its own kind than with the solvent 
(and more so with increasing discrepancies be- 
tween the molecular weights of rubber and 
solvent) ; and therefore the concentration in the 
rubber phase, even of the lowest molecular weight 
fractions, will always be higher than that in the 
solution phase, which is essentially solvent. 

As indicated above, this is a purely thermo- 
dynamic effect, and is not caused by mechanical 
co-precipitation or occlusion as has been sug- 
gested in the past. 

With increasing values of y», A increases,'* so 
that the tendency for the rubber to be in the 
precipitate phase increases; hence one obtains 
fractionation by systematically varying » in one 
direction or the other, taking advantage of the 
fact that the distribution between the two phases 
is strongly dependent on the molecular weight. 





=Bce. (43) 





(3) Comparison of the Precipitation and 
Extraction Methods 


It is now desirable to consider more in detail 
the two different methods of fractionation. ; 

In the precipitation method, the polymer is 
dissolved in a good solvent and then small 
increments of a non-solvent or precipitant (such 
as methyl alcohol) added until the solution 
breaks into two phases, with the formation of a 
swollen rubber phase or precipitate. This, of 
course, cannot occur until sufficient of the non- 
solvent is added to increase um beyond the critical 
value defined by Eq. (6). 

The original rubber of distribution function 
f(m) has now distributed itself between the two 
phases, and the distributions in the two phases 
are defined by f,(m) and f.(m) in Eqs. (37) and 
(38). The solution is then decanted and fraction- 
ation continued on it by adding more non- 
solvent and so increasing yw» further. Two phases 
exist again, with similar equations holding: 











i f.(m) 
iat 1+ ’c exp (A’m) 
f(m) 
- ; (44) 
(1+ ’c’e4’™) (1+ 8ce4™) 
ay fu(m) 
asi 1+[exp (—A’m) ]/p'c 
f(m) 
45) 





~ {1-+[Lexp (—A’m)1/B’c} (1-+Bce4™} 


More generally one may express the distri- 
bution in the mth fraction as: 





f."(m) = 





fp"(m) = 


The distribution function of the mth fraction is 
therefore represented by f,"(m). Actually, this 
may be very nearly represented by neglecting all 
but the last two terms in the denominator. 














f(m) 
46 
(1+Bac exp (Aam))(1+sc exp (Asm))- ++ (1+Bnc exp (Aqm)) (46) 
f(m) 3 si 
exp (—A,m) 
(1++6xc exp (Aam))(1-+Bsc exp (Aum))-+- (14 , ) 
n€ 
fp"(m)= f(m) 
exp (—A,m) 
(14+Bnc exp (4mm) ( 1+ 
Poe 48) 


16 8 decreases with increasing yu, but this is more than 
offset by the increase in A. 


This function starts out at a very small value 
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for m=0, increases to a maximum and then 
decreases to zero for m= ~. 


1 





lim 


=" (14+Becexp (Anm)( 1+ 





exp =") 
Bre 
Bre 
~ (14+Bnc)(14+B.¢) 
1 


lim == (). 


exp (dial i 
(1+6n¢€ exp (Anm))( 1+ 





=B nC; 


(49) 





mM 





7 
(S50) 

In the extraction method, the polymer is ex- 
tracted with a blend of solvent and non-solvent, 
such that a small fraction of the total polymer— 
largely low molecular weight material—is re- 
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moved. Here again the distribution between the 
solution and the swollen rubber is defined by 
Eqs. (37) and (38). The solution is decanted, and 
the precipitate is further fractionated by ex- 
traction with a solvent mixture of lower pm. 

The following two items are the two essential 
differences in the two methods: (1) In the pre- 
cipitation method, yu is systematically increased ; 
in the extraction method, y» is systematically 
decreased. (2) In the precipitation method, the 
fraction removed is in the precipitate, and the 
solution is further fractionated while in the ex- 
traction method, the fraction removed is in 
the solution, and the precipitate is further 
fractionated. 

Following the same line of reasoning as in 
Eqs. (44) to (47), except in the reverse direction, 
one finds that for the mth fraction, in the case of 
the extraction method 


f(m) 





f."(m) = 


exp (—Aym) 


; (51) 








exp (—A ym) 
Coane raed (a 
Bc 


): -+(1+Bn,c exp (A»m)) 
Bye 


f(m) 





fp™(m) = 








(14 ed Co a) 7s (14 





exp (—A,»m) ) 
Bm 


The distribution function for this fraction is therefore represented by f,”(m) which may be closely 


represented by 


f(m) 





f"(m)= 


(1+8,¢ exp (Anm))( 14 


It is to be seen that to this degree of approxi- 
mation, this fraction of the extraction method has 
the same distribution as the mth fraction in the 
precipitation method, since Eq. (53) is the same 
as Eq. (48). 

We may therefore say that, to the first ap- 
proximation, the results of the two methods are 
identical, and that essentially the same fractions 
are obtained in both ways. 

This is not exactly true, however, for the exact 
Eqs. (47) and (51) are not identical. In the first 
place, A and 8 are actually dependent on the 
actual distributions in the mixture, which are 


: (53) 





exp (—A =) 


Bre 





quite different for the two methods. Therefore 
(A os eistecetion ca (A a) precipitation } (54) 
(Bm) extraction * (Bm) precipitation: (55) 


This of course will introduce differences, but 
this effect is distinctly secondary compared with 
a second factor which we may call the tail effect. 

If we take the limiting cases of Eqs. (47) and 
(51) for low molecular weights, we find that 


Precipitation: 
Bre 


(1+ Bc): --(1+8,c) 





lim fp"(m) = 


m->0 


(56) 
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Extraction: 


(Bac) ++ «(B,c) 
li Pa => . 
aia OO TOE 


Hence, we see that since Bc<1 the most sig- 
nificant difference is that in the extraction 
method, far less low molecular weight material 
will be found in the higher fractions than is the 
case for the precipitation method. 

It would seem therefore that, normally, the 
extraction method of fractionation should be 
preferable to the precipitation method. Certain 
objections to such a choice might be raised: 

(1) Because of the large and compact nature of 
the rubber phase, equilibrium is reached much 
more slowly in the extraction method. 

(2) If insoluble gel is present, its great swelling 
will cause a tremendous overlapping in the 
higher fractions. 

(3) Because of the great sensitivity of the 
solubility to slight changes in y, it is very difficult 
to choose solvent mixtures properly in the ex- 
traction method. 

It remains for the individual investigator to 
choose between the methods. In the opinion of 
the author, the precipitation method, if modified 
to the extent of washing the precipitate or 
reprecipitating the fraction (in order to remove 
the ‘‘tail’’), is adequate for all ordinary purposes. 





(57) 


(4) Effect of Temperature 


The temperature has no effect upon the general 
considerations of fractionation in the preceding 
pages. However, since » is known both theo- 
retically and experimentally, to be a function of 
temperature, decreasing as T increases,!’ it 
follows that the temperature is involved in any 
fractionation in practice. In general, during the 
process of fractionation, by either method, the 
temperature should be kept relatively constant or 
anomalous results may be found. 

A variation of the precipitation method is that 
of inducing precipitation by cooling. Since yp 
increases with decreasing temperatures, this is 
thoroughly consistent with the theory developed 
above. 


17 See Eq. (5), and for further discussion of the depend- 
ence of 4 on various factors, see M. L. Huggins, Ann. Acad. 
Sci. 44, 431 (1943). 
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Fic. 4. The effect of concentration on the sharpness of 
concentration. 


(5) Effect of Concentration 


The influence of concentration is illustrated 
most effectively by the approximate Eqs. (48) 
(53). If this function is plotted against c, the 
“‘concentration’’ of the fraction!® 


c=v,'/Vo, (58) 


the curves of Fig. 4 result. The curves represent 
the percentage of material of molecular weight m 
in a fraction taken between um and (um+.01), for 
four different concentrations, 10-?, 10-*, 10-4, and 
10-5, the values of » chosen so that the maximum 
for each concentration lies in approximately the 
same range. 


18 This term “‘concentration’”’ is a trifle unclear. It is the 
volume of the rubber in the precipitate as compared with the 
total solution. If, for example, one is fractionally pre- 
cipitating a 1 percent solution of polymer into ten equal 
fractions, each fraction is 0.1 percent and c=0.001. In the 
case of the extraction method, ¢ varies in the course of the 
extraction, since, at the beginning of fractionation, the 
precipitate phase contains the whole polymer. 
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Figure 4 shows that there are two factors in- 
volved in varying c. 

(a) For smaller values of c, the fractionation 
is sharper; there is less ‘‘tail” effect, and the 
maximum is higher.!® Had a fraction in the lower 
molecular weight range been chosen, the curve 
would have been still sharper. 

(b) With higher dilutions, that is, smaller 
values of c, the value of u necessary for precipita- 
tion of a given fraction is greater. For example, a 
dilution of 10-fold requires an increase of u, (and 
um) by 0.01 unit. In the case of the lowest 
fractions, this shift is more appreciable. 

From the foregoing discussion it is clear that 
from a theoretical point of view, c should be kept 
as small as possible to obtain the sharpest 
fractionation. Practically, however, one cannot 
carry such a solution too far. Not only does this 
entail working with great quantities of solution, 
but, because of the shift of «1 with concentration, 
it may be next to impossible to separate the 
lowest fractions without adding tremendous 
quantities of precipitant.”° 


19 The fact that the maxima in Fig. 4 correspond to only 
about 60 percent of the material of m = 1700 (MW = 170,000) 
is not significant. Had a wider range been taken (i.e., 
Hn—pm=0.02), such that the fraction was larger, the 
maxima would have approached 100. 

20Tt would seem that, since the lowest fractions are 
fairly sharp, even when precipitated from a concentrated 
solution, the best procedure would be to start fractionation 
in a very dilute solution, but after separation of the higher 
fractions, toconcentrate the material to facilitate separation 
of the lower fractions. This is, however, a complicated 
procedure. 
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In practice, one must strike a compromise. 
Most fractionations are carried out from solu- 
tions where the initial concentration of the total 
polymer is about 1 percent. For a fractionation 
into ten equal fractions, this corresponds roughly 
to c=10-*. This seems to give fairly good results, 
judging both from the experimental results, and 
from the c=10-* curve in Fig. 4. 


(6) Effectiveness of Fractionation 


Some idea of the effectiveness of fractionation 
can be seen from Fig. 4. The distribution of each 
fraction is fairly broad especially for fractionation 
carried out at high concentrations. Especially 
objectionable is the “‘tail” effect, for the presence 
of even a small amount of low molecular weight 
material can alter the number average molecular 
weight appreciably. 

A better fractionation can be obtained either 
by refractionating or by washing the precipitate 
with a fresh sample of the precipitating mixture ; 
either procedure serves to remove the lowest sub- 
fractions, while not materially affecting the rest 
of the fraction. This solution may then be added 
to the decanted supernatant solution and frac- 
tionation continued. Experiments show that 
these methods are reasonably effective in elimi- 
nating the ‘“‘tail.”’ 

The author wishes to thank Professor Hugh S. 
Taylor and Dr. Michael Magat for constant 
encouragement and many helpful suggestions, 
and the Rubber Reserve Company for releasing 
the material for publication. 
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The Influence of Molecular Flexibility on the Intrinsic Viscosity, 
Sedimentation, and Diffusion of High Polymers* 
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Expressions for the mean square separation of chain ends and modifications of the formula 
for an ideal coil are discussed. On the basis of these and of the hydrodynamic theory of intrinsic 
viscosity, an interpretation of the modified Staudinger rule is offered. It relates the exponent a 
of the molecular weight to a flexibility parameter p of the chain in a given solvent, varying 
between zero and one (Eq. (4), (5)). Recent data on polystyrene and on cellulose nitrate are 
analyzed in greater detail. By means of the frictional ratio f/fo, the sedimentation constant s 
and the diffusion constant D, respectively, are connected with the degree of polymerization in 
terms of p (Eq. (9)). The limiting dependence of sedimentation and diffusion rate upon mo- 
lecular weight for a straight chain and an ideal coil is also found in this manner. A comparison 
shows satisfactory agreement between values for p found from intrinsic viscosity and those 
determined from sedimentation or diffusion rates, for certain cellulose esters and starch deriva- 
tives. Effects of solvent and of inhomogeneity in respect to molecular weight are discussed 
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INTRODUCTION 


HE physical methods employed to deter- 
mine molecular weights of high polymers 
in solution are based either on the measurement 
of a colligative property like osmotic pressure or 
sedimentation equilibrium, or arate phenomenon, 
such as the viscosity increment due to the solute 
or the rates of sedimentation and diffusion of the 
solute. The second type of procedure is more 
easily applicable from the experimental point of 
view, while its interpretation is more involved 
than that of measurements of the former type. 
The determination of polymerization degrees 
from viscosity data is based on Staudinger’s 
empirical rule or a modification thereof obtained 
by calibration of the procedure by means of an 
absolute method. In this way an empirical 
relation is obtained between the intrinsic vis- 
cosity and the molecular weight. In the case of 
sedimentation-diffusion, an expression for the 
molecular weight is derived in a straightforward 
way. For polydisperse systems however, it is not 
immediately evident what type of molecular 
weight average results. 
Present theories on the other hand, establish 
relations between the intrinsic viscosity and the 
shape of the solute molecule.! Similarly, the result 


* Publication assisted by The Graduate School of Howard 
University. ; 
1 See for instance R. Simha, J. App. Phys. 13, 147 (1942). 


of a sedimentation-diffusion measurement can be 
related to the shape of the sedimenting particle 
by means of the so-called frictional ratio.* A con- 
necting link between some of the empirical 
results in terms of molecular weight and present 
theories can be established, if the shape of long 
chain molecules in solution can be related to 
their degree of polymerization. In this manner 
it becomes possible also to relate the intrinsic 
viscosity and the molecular weight to the indi- 
vidual sedimentation or diffusion constants in- 
stead of their ratios only. It is the purpose of 
this article to examine certain aspects of these 
questions. 


I. SHAPE OF LONG CHAIN MOLECULES 
IN SOLUTION 


It has been shown by various authors* follow- 
ing Rayleigh’s‘ treatment of the problem of 
random flight, that the mean square separation 
of ends in a chain consisting of N freely rotating 
links is given by 


Jim (R?)w=AN, (1) 


where A is the squared length of a link and R the 


2 The Uliracentrifuge, Ed. The Svedberg and Kai O. 
Pedersen (Oxford University Press, New York, 1940). 

3W. Kuhn, Kolloid Zeits. 68, 2 (1934); E. Guth and 
H. Mark, Monats. 65, 93 (1934). 

4 Compare the excellent discussion by S..Chandrasekhar, 
Rev. Mod. Phys. 15, 2 (1943). 
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distance between chain ends. The constancy of 
the valence angle merely introduces a numerical 
factor. Relation (1) will remain valid also if an 
energy barrier to rotation exists, provided that 
rotational symmetry is preserved (ethane-type 
of barrier). The factor A is modified and be- 
comes a function of the height of the energy 
hump and of the temperature.’ Equation (1) 
can be generalized. If a distribution of ‘“‘step”’ 
lengths, 
r(x, y, 2)dxdydz, 


exists in proceeding from the ith joint to the 
(i+1)th, and 7 is independent of the ordinal 
number 7, where x, y, z denotes a set of coordi- 
nates through the joint, one obtains‘ 

lim (R?)y=AN+BN?. 


N-o 


(1a) 


A indicates the mean square of the length of a 
“step” or effective length of a link, averaged over 
the distribution 7. B denotes the sum of the 
squares of the mean displacements in 3 directions 


B=(EP+(nP +5 /. 

The system é, 7, ¢ results from the original one 
by means of a certain orthogonal transformation 
which need not be specified here. If the function 
t possesses the previously mentioned symmetry 
properties, B=0. Thus again a random distribu- 
tion of chain end distances is obtained, however 
superimposed upon that resulting from a straight 
extension of the chain. B is 4 measure of the 
geometrical anisotropy existing in the structure. 

Considering now an actual chain molecule, the 
fact that ‘‘tracing back of steps’’ or back coiling 
of the chain is restricted because of the space 
filling effect of the substituents, requires a modifi- 
cation of the considerations leading to rela- 
tion (1). It is extremely difficult to evaluate this 
effect in detail, since the direction of each suc- 
cessive link now becomes dependent upon that 
of the preceding one. However, it may be seen 
that, on the average, a straightening out of the 
chain coil will occur, larger values of the separa- 
tion of chain ends, than in the ideal case, being 
more probable. The actual chain may therefore 
be replaced in a first approximation by one in 
which there exists a distribution 7 of displace- 
ments around each joint, and independent of its 


5S. E. Bresler and J. I. Frenkel, Acta Physicochimica, 
U.S.S.R. 11, 485 (1939). 
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position within the molecule, giving rise to a 
mean value (R?*), as in (1a). For the large degrees 
of polymerization involved, even slight devia- 
tions from symmetry can lead to appreciable 
values for the second term in (1a). Unless the 
range of N-values is large, this equation can be 
approximately replaced by a single term with an 
intermediate exponent, thus 


(R®)w=CN't*#; 0<x<1, (1b) 


C is a constant, «=0 corresponds to the ideal 
coil, as defined in (1), x=1 to a straight chain. 
Thus we arrive at a result similar to the one 
suggested previously by Kuhn.’ It may be noted 
that in the transition from (1a) to (1b), x in- 
creases with B while C decreases. For N ~10', 
a value of x of 0.4-0.5 can increase the value of 
((R?))* for a paraffin chain by a factor of several 
units. Finally the influence of the solvent must 
be taken into account. Stuart’s* experiments on 
models propose to reproduce the difference be- 
tween an inert and a good solvent. In the latter 
system the chains are straightened out. The 
solvent in this instance prevents the polymer 
molecule from curling up. This may be due in 
part to actual solvation and in part may be 
considered analogous to the conditions which 
prevent the formation of a drop and lead to 
spreading of a film. As is well known,’ this 
occurs when the energy of adhesion between 
upper and lower liquid exceeds the cohesion 
within the former one. Of course, in our case, 
the entropy factor is the primary source for the 
tendency in the single polymer molecule to form 
a “drop.’”’ Furthermore these considerations can 
apply to dilute solutions only. 


II. INTRINSIC VISCOSITY, MOLECULAR 
WEIGHT, AND DIMENSIONS 


Since Staudinger first proposed his relation, it 
has become necessary to modify it. For relatively 
small molecular weights, not exceeding 30,000 as 
a weight average, addition of a constant term® 
has been found satisfactory. For larger molecular 


6H. A. Stuart, Naturwiss. 31, 123 (1943). 

7 Compare for instance N. K. Adam. The Chemistry and 
Physics of Surfaces (Oxford University Press, New York, 
1941). 

8 R. Fordyce and H. Hibbert, J. Am. Chem. Soc. 61, 
1912 (1939); P. J. Flory and P. B. Stickney, J. Am. Chem. 
Soc. 62, 3032 (1940); W. O. Baker, C. S. Fuller and J. H. 
Heiss, Jr., J. Am. Chem. Soc. 63, 3316 (1941). 











weights, in accordance with Houwink,? the fol- 
lowing equation has been applied 


ni=KM*=K’'Z:, (2) 


where 7; is the intrinsic viscosity, found by 


extrapolation from a plot of 9.)/c vs. c. M is the 
molecular weight of the fraction as determined by 
an absolute method, Z the degree of polymeriza- 
tion. Since the partial specific volume of the 
polymer is almost always practically independent 
of molecular weight, the result remains unaltered 
if c is expressed as a volume ratio rather than a 
ratio of weight to volume. A few examples are 
presented in Table I. The molecular weights were 
determined partly from osmotic pressure and 
partly from ultracentrifuge data. In several 
cases only three points were available for com- 
puting a. The homogeneity of the samples varied, 
those investigated by Mosimann and Gralén 
being quite inhomogenous. The meaning of the 
last column is discussed later on. This table is by 
no means complete. However, it seems to indicate 
that the exponent a lies roughly between 0.5 
and 1.5. It may be noted that relatively slight 
inaccuracies in the viscosities or’ molecular 
weights can produce noticeable variations in a. 

As mentioned before, present hydrodynamic 
theories of intrinsic viscosity introduce as funda- 
mental parameter the shape of the solute mole- 
cule.! If it is approximated by an ellipsoid of 
revolution, the length-thickness ratio //d appears 
in the final result.! For the case of overwhelming 
Brownian motion in the viscosimeter, »; becomes 
a rather complicated function of //d, which has 
been evaluated in the range between //d=1 and 
300." Using these data, Kraemer” has given an 
interpolation formula in the form of a polynomial. 
However, a simpler equation can be given, which 
is more useful for our purposes, namely : 


ni~0.233(1/d)8: 20<1/d<100; 

Average deviation 1.2 percent (3) 
ni ~0.207(1/d)'™; 2051/d=300; 

Average deviation 2.3 percent 


®R. Houwink, J. prakt. Chem. 157, 15 (1940). 

10W. Kuhn, Helv. Chim. Acta 26, 1394 (1943), where 
previous references are given. In the present author’s 
opinion, Kuhn’s hydrodynamic treatment and others de- 
rived therefrom, contains an incorrect assumption for the 
limit of vanishing shear gradient, which makes the result for 
ni too small. See R. Simha, J. Phys. Chem. 44, 25 (1940). 

J. W. Mehl, J. L. Oncley and R. Simha, Science 82, 
132 (1940). 
2 E. O, Kraemer, J. Frank. Inst. 231, 1 (1941), 
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SIMHA 
TABLE I, 

Solute Solvent a p 
Polyisobutylene Diisobutylene 0.642 0.38 
Polystyrene Toluene 0.70° 0.42 

(polymerization 

temperature 60°C) 
Polystyrene, 120°C Toluene 0.81° 0.47 
Polystyrene, 180°C Toluene 1.145 0.68 
Cellulose Schweitzer’s reagent 0.80° 0.47 
Amylose acetate Methyl acetate 0.954 0.57 
Cellulose acetate Acetone 1.0° 0.59 
Cellulose nitrate Acetone 1.0/ 0.60 
Polyvinyl chloride Cyclohexanone 1.09 0.59 
Amylose Ethylene diamine 1. 0.88 








@P, J. Flory, J. Am. Chem. Soc. 65, 372 (1943). 

+ See reference 15. 

¢ See reference 19. 

4 See reference 22. 

¢M. Harris and A. Sookne, private communication. Molecular 
weights from osmotic pressure. 

4 See reference 18. 

9 D. J. Mead and R. M. Fuoss, J. Am. Chem. Soc. 64, 277 (1942). 

4 J. F. Foster and R. M. Hixon, J. Am. Chem. Soc. 66, 557 (1944). 


With increasing axis ratio, the exponents increase 
slowly but do not reach a value of two within the 
range of importance. For smaller values of //d, 
the exponent in (3) assumes a value of 1.72, the 
factor equals 0.217, and a constant equal to 2.283 
is to be added, similarly as in the correction to 
Staudinger’s rule.* »; is expressed in terms of 
volume concentration. The axis ratios computed 
from measured viscosities represent the dimen- 
sions of the effective ellipsoid which makes the 
same contribution to 7; as the actual molecule. 
This is certainly not more than a rough approxi- 
mation. For high degrees of polymerization, 
where the chains form preferentially compact 
coils,” this should be more satisfactory than for 
intermediate ones. For short chains on the other 
hand, where curling up is insignificant, this 
approximation should again be more adequate. 
However in the following, we shall be concerned 
only with the former extreme. Furthermore, the 
possibility of solvation in good solvents con- 
tributing to the observed large intrinsic viscosities 
has been disregarded and the whole effect at- 
tributed to the asymmetry of the unsolvated 
structure. In the case of globular proteins, it has 
been possible to discriminate to some extent 
between these two possibilities by comparison of 
viscosity and diffusion-sedimentation data." The 
difficulties encountered at present in the evalua- 
tion of such measurements for chain polymer 
systems, make this possibility appear doubtful in 
our cases. However, for high asymmetries, the 
solvation effect can produce only errors which 
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Fic. 1. Double logarithmic plot of asymmetries //d vs. 
degree of polymerization Z for polystyrene: 


I Polymerization temperature 60°C. 
II Polymerization temperature 120°C. 
III Polymerization temperature 180°C. 


are minor compared with those arising from other 
difficulties. In our applications two further points 
must be kept in mind. First the fact, that for a 
given degree of polymerization, the length and 
diagonal dimensions are only statistically defined. 
Therefore Eq. (3) must be modified to contain 
((1/d)*)4, where q is the exponent in (3) and the 
average is formed over the distribution of dimen- 
sions in the coil in solution. The exact form of 
this distribution is not known except for ideal- 
ized cases, in which expressions of an essentially 
Gaussian type are obtained.‘ However, the above 
average will differ from ((//d)?)", only by a 
numerical factor of the order of unity. Secondly, 
the polydispersity of the solute must be taken 
into account. The mean axis ratio determined 
from the intrinsic viscosity of a mixture is: 


(1/d)m= {20 (1/d) Wi}, 


where W;, is the weight fraction of species 7 with 
axis ratio (//d); in the mixture. The higher the 
molecular asymmetry, the greater the effect of 
inhomogeneity.” 

To provide, with these reservations, a bridge 
between Eqs. (2) and (3), the dimensions of the 
polymer molecule must be expressed as a function 
of the degree of polymerization. It is not difficult 
to obtain the result in two limiting cases, namely 
those of a perfectly straight chain and an ideal 
flexible coil respectively. In the first instance, the 
longest dimension increases proportionally to 
the molecular weight, the cross section remains 
constant in a homologous series, and the axis 


ratio is therefore proportional to Z. In the other 
extreme, the mean separation of chain ends is 
given by relation (1). Kuhn* has shown that in 
such a structure on the average all dimensions 
increase at the same rate with the molecular 
weight. Thus the axis ratio becomes independent 
of the degree of polymerization. Our previous 
considerations in Part I lead for a chain in solu- 
tion to an equation of type (1b), corresponding 
to a certain anisotropy of the structure. We may 
therefore expect an expression,! 


l/d=CZ?; 0<p<i, (4) 


to hold approximately. The lower limit of the 
exponent characterizes the ideal coil, the upper 
one a straight rigid chain. Combination of Eqs. 
(4) and (3) leads directly to the empirical for- 
mula (2) with 

a+=1.7p. (5) 


In this manner most of the values in the last 
column of Table I were computed, with the 
exception of the values for polystyrene and 
nitrocellulose which were found directly from 
suitable plots and are slightly different. However. 
(5) is sufficiently accurate for a rapid calculation 
of p from a or vice versa. It follows from (4) and 
(5), that the maximum. possible value of the 
viscosity exponent a should lie around 1.7. In 
those cases in which Staudinger’s equation holds, 
a value for p of 0.59 results, approximately half- 
way between zero and one. According to this 
theory no particular physical significance is to be 
attached to the case a=1. Furthermore it is 
understandable for instance, that for polyiso- 
butylene the value lies below that for the cellulose 
esters. Also it may be expected that for a given 
polymer, in a poorer solvent or solvent mixture, 
the p-values will be smaller than in a thermo- 
dynamically better one. The behavior of cellulose 
acetate in acetone-methanol mixtures seems to be 
in agreement with this conclusion, since p de- 
creases with the addition of methanol.’* Measure- 
ments of the depolarization factor of the light 
scattered by the particle'* should also reveal an 
increased degree of coiling upon addition of 
precipitant. Furthermore, preferential forward 
scattering should appear at higher molecular 
weights in the poorer solvent system. 


13H. Mark, private communication. 
4 See for instance P, Debye, J. App. Phys. 15, 338 (1944). 
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Fic. 2. Double logarithmic plot of asymmetries //d vs. molecular 
weight M for cellulose nitrate. 


Figure 1 gives logarithmic plots of //d, as de- 
termined from viscosity,!*" vs. Z for the three 
series of polystyrenes. Equation (4) is satisfied. 
Tables II-IV contain the experimental data to- 
gether with the dimensions / and d computed 
therefrom. No data on the partial specific vol- 
umes of the fractions were reported. Therefore a 
value of 0.93 was assumed for all three series. 
The third column was computed from Eq. (4). 
Using a value of 2.5A-units for the monomer, 
a maximum length /n.x can be found. Jmax/l varies 
in the three series respectively from 15.5 to 9.1, 
6.3 to 4.2, and 5.4 to 4.9, indicating degrees of 
coiling increasing with Z, as is to be expected. 
An increase of the viscosity exponent a and the 
flexibility parameter p is accompanied by an in- 
crease of the slope of the reduced osmotic pressure 
curves.'5 According to present theories,'® this 
indicates a better solvent, in which a relative 
straightening out of the chains becomes under- 
standable. On the basis of the fact, that his 
‘constants’ K decreased with increasing tem- 
perature of polymerization, Staudinger suggested 
that branching is favored under such conditions. 
This question can of course not be decided on the 
basis of viscosity alone. True, the proportionality 
factors in the equations for //d decrease between 

16 T, Alfrey, A. Bartovics, and H. Mark, J. Am. Chem. 
Soc. 65, 2319 (1943). I am indebted to Dr. H. Mark for 
making available the data in more detail than in the 
publication. 


1 P, J. Flory, J. Chem. Phys. 10, 51 (1942); M. L. 
Huggins, Ann. N. Y. Acad. Sci. 43, 1 (1942). 


60° and 180°. However, this may not be more 
than a consequence of the increase in p, as 
discussed in connection with Eqs. (1a) and (1b). 
This increase in p indicates increased stiffening 
of the chains. Now a branched structure offers 
more resistance to convolutions of the chain than 
a linear one. Furthermore the chains become 
more symmetrical (smaller //d), for the same 
degree of polymerization, as may be seen by 
comparison of the first and last rows in Table 1V 
with the last two in the preceding one. Clearly 
more measurements on overlapping fractions 
prepared at different temperatures are needed.!” 

Mosimann’s!* measurements on nitrocellulose 
were also analyzed. Figure 2 shows plots of 
log 1/d vs. log M. Equation (4) is applicable to the 
four highest molecular weights, ranging from 
613,000 to 30,000. The degrees of nitration are 
slightly different. The two smallest degrees of 
polymerization were found by means of equi- 
librium measurements in the ultracentrifuge, the 


TABLE II. Polystyrene 60°. 











Z(obs) Z(calc) 1x10-3 d 

ni x10-3 x<10-3 l/d (A) (A) 
355 19.7 19.6 75 3.18 42 
295 15.2 15.0 67 2.77 41 
200 8.6 8.8 53.5 1.95 36 
140 5.3 5.2 43 1.45 34 


ni =3.37XK1072Z-7; 1/d=1.17Z2-"; 1=7.50Z-°%; d=6.41Z-9 








( 17 See also M. L. Huggins, J. Am. Chem. Soc. 66, 1991 
1944), 
18H. Mosimann, Helv. Chim. Acta 26, 61 (1943). 

















‘in the manner described below. 
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remainder from rate data. No inhomogeneity 
coefficients are given for these samples except for 
M=30,000 which had a ratio of 1.65 between 
weight- and number-average molecular weight. 
In view of the magnitude of the average axis 
ratios, ranging from 285 to 47.5, this represents 
quite high a value. In the original article, values 
for the molecular asymmetries are presented. 
They were computed on the basis of theories 
criticized elsewhere” by the present author. 
A comparison of these //d-values with the maxi- 
mum values for a straight chain also shows 
that they are too large. Using our values for //d, 
ratios of Jmax/l between 2.7 and 1.2 were com- 
puted. No detailed analysis of other data pre- 
sented in Table I was undertaken but merely 
approximate values of » were computed from 
Eq. (5). 

Equation (4) is based on (1b). For wider 
ranges of molecular weights in a series, it may 
be necessary to use a formula modified in accord- 
ance with (la). The Staudinger-type relation 
(2a) would then be replaced by a corresponding 
two-term equation. 


Ill. SEDIMENTATION, DIFFUSION, AND 
MOLECULAR WEIGHT 


The basis for the determination of molecular 
weights from rate measurements in the ultra- 
centrifuge is given by the familiar equation?: 


M=(RT/1—Vop)-s/D. (6) 


(1— Vp) represents the Archimedes-factor, V in- 
dicating the partial specific volume of the solute, 
p the density of the solvent. s is the sedimenta- 
tion constant, D the diffusion constant of the 
particle. The importance of extrapolating both 
parameters to zero concentration has been re- 
peatedly emphasized.'!*!* These constants can be 


TABLE III. Polystyrene 120°. 











Z(obs) Z(calc) 1x<10-3 d 
ni x<10-3 x<10-3 l/d (A) (A) 
155 3.3 HK 46.3 1.29 28 
118 as 2.3 39.5 1.03 26 
86 1.6 1.6 33.0 82 25 
62 1.0 1.1 27.0 .62 23 


n=2.26X 1072; 1/d=1.02Z 47; 1=6.83Z-; d=6.72Z 18 








19N. Gralén, Thesis, Upsala, 1944. The exponent is 
calculated from the author’s sedimentation-diffusion data 





related to the shape of the sedimenting particle 
by means of the frictional ratio f/fo.? It indicates 
the ratio between the frictional resistance of the 
actual molecule and a sphere occupying the same 
volume. If the partial specific volume V of the 
polymer is independent of the molecular weight, 
the following equation is obtained : 


f/fo=const./DM}. (7) 
For an ellipsoidal particle, f/fp can be expressed 


TABLE IV. Polystyrene 180°. 








Z(obs) Z(calc) 1x10-3 d 





ni x10-3 x10 I/d (A) (A) 
73 1.6 1.6 29.5 75 26 
59 1.4 1.4 26.0 .67 26 
53 1.3 1.3 24.5 62 25 
45 1.1 1.1 22.0 54 25 
ni=1.58X1077Z!-4; 1/d=1.91K 102-8; 


1=2.252-%; d=11.74Z" 








as function of //d.?° It has been tabulated up to a 
value of hundred for the axis ratio,? and recently 
for higher values up to three hundred.” Although 
the original equation is not as cumbersome to 
evaluate as the viscosity formula, it is again 
useful for our purposes, to write it as an ex- 
ponential expression. The interpolation formula: 


f/fo= 0.506(1/d)*%; 20S1/d=300; (8) 
Av. dev. 1.3 percent, 


represents the values well. 

Axis ratios calculated from the frictional ratios 
and those obtained from viscosity measurements 
have been found to agree!'*! in the case of 
proteins. The discrepancies become greater for 
chain polymers, as indicated, for instance, by 
measurements of Signer and collaborators,”! 
Gralén,’ Dombrow and Beckmann” and of 
Mosimann.!* A comparison of Eqs. (8) and (3) 
shows that for polydisperse systems, different 
types of averages for //d are obtained. Further- 
more the molecular weights in (7), if determined 
from (6), are usually calculated as a ratio of 
some average sedimentation constant to a weight 
average diffusion constant, rather than as an 
average ratio of the two constants. Finally the 


20F. Perrin, J. de phys. et rad. 7, 1 (1936). 

1M, A. Lauffer, J. Am. Chem. Soc. 66, 1188 (1944). 

2B. A. Dombrow and C. O. Beckmann, private com- 
ae, Molecular weights from sedimentation-dif- 
usion. 











194 ROBERT SIMHA 

















7 
G 
s 

4r 
Q3- 

e 
o 

er 

4 r 4 i 
ts 30 sO 


20 
49-5 


Fic. 3. Double logarithmic plot of diffusion constant D 
vs. molecular weight M for methyl cellulose. Circles indi- 
cate experimental values, dashed line the experimental 
slope, solid line the computed one. 


accuracy in determining f/f) is limited to a 
greater extent than that in viscosity measure- 
ments, errors up to 20 percent not being out of 
the ordinary. This makes the above-mentioned 
discrepancies understandable. 

It is assumed now that the particle in the 
ultracentrifuge has the same //d-ratio as in the 
viscosimeter for a given solvent. This should be 
true for small shear gradients, since the deforma- 
tion of a flexible coil is a second-order effect.’ 
Then it is possible to derive an expression for the 
separate dependence of the sedimentation or 
diffusion constant upon the molecular weight 
from viscosity data by elimination of the axis 
ratio //d. Introduction of formula (4) into (8) 
and combination with (7) yields the result: 


D=const M-‘#-4547).  s=const Mt-454?, 

Consequently for p=1: 0) 
D=const M79, s=const M°*!; 

for p=0: D=const M-?, s=const M}. 


The constants are independent of the molecular 
weight in a given series, having a given p-value. 
In the limit of a straight chain, the sedimenta- 
tion constant becomes almost independent of 
molecular weight. This is to be expected, since s 
depends for long chains mainly upon the ratio 
M/l. If Staudinger’s rule holds, p=0.59 and 


D=const M~-®, s=const M*”, (9a) 


Equation (9) is not strictly applicable to the 
extreme of an ideal coil (p=0), since it is based 
on relation (7), valid if the molecular volume of 
the polymer increases with the first power of. M. 


For the average Kuhn-coil, however, it is pro- 
portional to M}, thus causing the diffusion con- 
stant to vary with the reciprocal square root 
of M. The experimental data presented below 
give exponents in the neighborhood of 3. Their 
interpretation on such a basis, however, is in- 
compatible with the observed viscosity behavior. 
Kuhn’s recent considerations” are not in agree- 
ment with (9) and (9a). This arises from the 
corresponding divergence in the interpretation 
of Staudinger’s rule, pointed out in Part II. 
A formula of type (9) has been proposed by 
Gralén,!® without any relation to viscosity or the 
degree of internal flexibility of the chain. From 
his empirical exponents for cellulose, cellulose 
nitrate in acetone and sodium cellulose xanthate, 
approximate values of p equal to 0.47, 0.78, 0.41 
and of a equal to 0.80, 1.32 and 0.70 are cal- 
culated. 

For an inhomogenous system, Eq. (9) indi- 
cates that no number- or weight-average mo- 
lecular weight is obtained by forming the ratio 


(s)/{D)w 


The same holds for measurements of free diffu- 
sion or of ultracentrifuge rates. Equation (9) 
indicates what type of average is to be expected, 
if the viscosity-molecular weight relation (2) is 
known for the series. It may be noted that the 
variations in the dependence of s and D upon M 
with the nature of the polymer and solvent 
system are much slighter than in respect to 
viscosity behavior. 

In Table V are presented p-values computed 
from the observed dependence of s or D upon M 
(Eq. 9), and are compared with those obtained 
from viscosity by means of relation (5) or di- 
rectly from the axis ratios. In Figs. 3 and 4, 


TABLE V. Values of flexibility parameter p. 











From Max. 
From sedimen- From dev. 
Compound diffusion tation viscosity % 
Methy] cellulose 0.57 — 0.59¢ 5 
Cellulose nitrate? — 4 0.60 11 
Amylose acetate? 0.48 0.49 0.53-0.57/ 11 (sed.) 








@ See reference 2. 

+ See reference 24. 

¢ The values for //d originally given in Table III of reference 1 are 
erroneous and should be replaced respectively by 38.5, 55.5, 69.5, 
resulting in the above exponent. 

4 See reference f of Table I. 

¢ See reference 22. 

J The lower value results from a comparison of the axis ratios deter- 
mined by means of f/fo with the molecular weights. The upper one is 
derived from the viscosity values of //d and has been used in Fig. 4. 
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Fic. 4. Double logarithmic plot of sedimentation constants s vs. molecular 
weight M for amylose acetate (A) and cellulose nitrate (C). Circles indicate 
experimental values, dashed lines experimental slopes, solid lines com- 


puted ones. 


D or s are plotted vs. M for these three polymer 
systems. For cellulose nitrate, the slope of the 
line through the second and the last point was 
calculated ; in the other cases the two end points 
were used. The computed D or s were fitted 
through the second point in the former case 
and through the first one otherwise, using the 
p-values from viscosity in Eq. (9). The extent of 
the agreement between viscosity data and ultra- 
centrifuge results, as interpreted by the present 
theory, is shown by a comparison of the slope of 
the solid lines (the » from viscosity data), and 
that of the dashed lines (the » from rate data). 
The last column in the table above contains the 
maximum deviations between experimental sedi- 
mentation or diffusion data and those calculated. 
In view of the accuracy obtainable and of factors 
discussed above, it does not appear necessary to 
determine the slopes in a more accurate manner. 
The molecular weights used in Figs. 3 and 4 
represent ultracentrifuge values. It would be 
very desirable to have osmotic determinations in 
addition. The agreement obtained seems satis- 
factory, in spite of the fact that in two cases only 
two points were available for determining the 
slope, as may be judged from the graphs. At- 
tempts to apply this theory to Signer’s and 
Gross’** sedimentation data on some polystyrene 
samples and to Polson’s** data on cellulose 

23 R. Signer and H.Gross, Helv. Chim. Acta 17, 348 (1934). 


24 Computed from data by A. Polson, Kolloid Zeits. 83, 
172 (1938). 


acetate meet with failure, the exponents p turn- 
ing out to be too small or too large. .It would be 
very desirable to have sedimentation and diffu- 
sion data on the polystyrene and cellulose acetate 
fractions referred to in Table I. 

Summarizing, it may be stated that by means 
of a combination of the present theory"! of 
intrinsic viscosity of long molecules with certain 
semi-empirical results on the shape of flexible 
chains in solution, it is possible to interpret the 
empirical relation between intrinsic viscosity and 
average degree of polymerization arising from a 
modification of Staudinger’s rule, and to account 
for variations in the exponent with the nature of 
the polymer and of the solvent or solvent mix- 
ture.” Application of the hydrodynamic theory of 
sedimentation and diffusion combined with the 
above results gives the sedimentation or diffu- 
sion constant as a function of the molecular 
weight and the variation of this function with 
the internal flexibility of the chain and with the 
medium. 

The author wishes to express his appreciation 
to Dr. John J. Grebe for the interest shown in 
this work. 

Aid in carrying out the numerical calculations 
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% The influence of such factors on higher terms in the 
viscosity-concentration relation has since been discussed 
by E. M. Frith, Trans. Faraday Soc. 41, 17, 90 (1945). 
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The mechanical laws of impact lead to the conclusion that, in any transfer of kinetic energy 
into quantized energy (excitation, dissociation), the “excitation function’ should, at the 
threshold, start from zero and increase with increasing kinetic energy. This shape of the 
function is likely to cause a systematic error in the measurement of critical energies with 


certain non-equilibrium methods, 





I, PROBLEM 


HE procedures for the determination of 
chemical energies, for example, energies of 
dissociation, may be classified into equilibrium 
and non-equilibrium methods. Chemical egui- 
libria are not affected by probabilities of transi- 
tion between the various states present, since, in 
the underlying thermodynamic relations, only 
properties of states are involved (energies, sta- 
tistical weights), but no probabilities of transition. 
The situation is different for non-equilibrium 
methods in which the rate of a process is ob- 
served as a function of the temperature. The 
accuracy of such methods is restricted by lack 
of information regarding probabilities of transi- 
tion. This will be evident from an example 
recently discussed by Dwyer and Oldenberg! in 
their comparison of the various methods for the 
determination of the energy HXO-H+OH. We 
consider the breaking up of H2O molecules into 
H+OH by impacts of the second kind with 
excited mercury atoms. Since the energy of 
excitation of the mercury atom is slightly de- 
ficient, the energy of the thermal collisions 
provides the balance if the temperature is high 
enough. Thus the increase of rate observed with 
increasing temperature is attributed to the in- 
crease of the number of thermal collisions occurring 
with sufficient energy. In this argument the same 
effect is attributed to any thermal collision 
occurring with an energy above the threshold. 
Hence the argument fails to take into account 
the possibility that a collision occurring with an 
energy slightly above the threshold may have an 
efficiency strongly dependent upon this excess 
1R. J. Dwyer and O. Oldenberg, J. Chem. Phys. 12, 353 


(1944). See O. Riechemeier, H. Senftleben, and H. Pastorff, 
Ann. d. Physik 19, 202 (1934). 


energy. This dependence (efficiency of collisions 
against excess energy) is called the excitation 
function of the process. The excitation function 
must be known for a quantitative, theoretical 
interpretation of the rate observed as a function 
of the temperature. 

The only process for which excitation functions 
have been measured is the excitation of atoms by 
electron impact. In this instance the efficiency 
of the impact, starting from zero at the threshold, 
rapidly increases with the excess kinetic energy 
of the impinging electron.’ 

In order to judge non-equilibrium methods we 
must find out whether the rapid increase of the 
excitation function just above the threshold is a 
special property of the process just mentioned or 
is expected to apply as well to other processes. 
A special case which can be analyzed in some 
detail leads to a general conclusion regarding 
excitation functions right above the threshold. 


Ii. SPECIAL CASE 


As a special case we analyze the excitation of 
molecular vibration by electron impact. Accord- 
ing to Eucken and Franck,’ the molecule simply 
picks up vibrational energy, while only very 
small amounts of linear or angular momentum 
can be transferred because of the ratio of masses 
of electron and molecule. 

The “controlled electron impact experiment” 
deals with impacts in which a fast electron hits 
an atom or molecule which is at rest with respect 
to the observer. Such impacts fail to cause 
quantum transitions if the kinetic energy W has 


2 W. Hanle and K. Larché, Ergeb. d. exakt. Naturwiss. 
10, 285 (1931). 

3 J. Franck and A. Eucken, Zeits. f. physik. Chemie B20, 
460 (1933). 
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merely the value Q of the energy of excitation 
computed from spectroscopic data. This failure 
is due to the fact that, according to the laws of 
conservation, the common center of gravity of 
the system continues its path after the impact so 
that the kinetic energy of the impinging particle 
is not fully available for the quantized process. 
In the impact most favorable for excitation, that 
is the central impact, W must at least equal 
Q(M-+m)/M. (M and m=masses of atom and 
impinging particle, respectivély.)* 

For excitation by impact of electrons the factor 
(M+m)/M is nearly unity, so that no noticeable 
excess of the kinetic energy W over the excita- 
tion energy Q is needed for the central impact. 
We are interested, however, in another aspect of 
the impact. For off-center impacts the electron 
will bounce off obliquely and retain a larger part 
of its initial energy. Therefore, in this case, only 
an electron with higher kinetic energy W will 
be able to produce excitation. Vice versa, a 


faster electron has a greater chance for being 
effective because it can act anywhere within a 
larger central part of the cross section of the 


target atom. Thus we can correlate the kinetic 
energy of the impinging electron with the proba- 
bility of excitation. That is, we can make a 
prediction regarding the excitation function. 
This function must start from zero because, for 
the electrons which have exactly the threshold 
energy, the effective cross section is vanishing.® 

We do not claim that this kinetic theory 
furnishes a complete explanation of the excita- 
tion function. Instead, the complete representa- 
tion will consist of two factors, the kinetic factor 
discussed here and a structure factor, which is not 
accessible to a mechanical treatment. 

Let us consider an electron incident on the 
spherical surface of the molecule with the kinetic 
energy W under an angle a (Fig. 1). It bounces 
off under another angle, delivering the quantized 
energy Q. The computation of this process is 
based on the laws of conservation. The equations 


* On the basis of this restriction G. Joos and H. Kulen- 
kampff computed the minimum kinetic energy of positive 
ions required for the excitation of atoms by impact 
(Physik. Zeits. 25, 257 (1924)). Qualitatively, this restric- 
(192 DP discussed by J. Franck (Zeits. f. Physik 25, 312 

5 This follows from the equations from which Klein and 
Rosseland concluded the occurrence of impacts of the 
second kind (Zeits. f. Physik 4, 49 (1929)). 
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so obtained, which may here be omitted, have 
real solutions only if cos? a=Q/W, that is, if the 
kinetic energy of the electron associated with the 
normal component of the velocity is equal to or 
larger than Q. Thus, for a given kinetic energy 
W of the electron, it is possible to compute a 
limiting value of the angle of incidence @ and, 
hence, the central section mr? (see Fig. 1) of the 
target molecule within which the electron of the 
energy W possibly is effective. This section 
divided by the total area 7R? of the target 
represents the effective fraction of the target, 
that is, the kinetic factor k of the excitation 
function. The result is 


k=1-Q/W [for W2Q]. 


As expected, this simple function represents a 
curve starting from zero at W=Q, approaching 
unity for large W. 

The same computation carried through for the 
case of comparable masses m and M leads to 
the more general equation 


QO 4M+m M+m 
b=1——( fo w=0 —*)|. 
W\ M M 


In these computations the molecule has been 
represented as a sphere with a hard surface. But 
the theory should apply as well to a molecule the 
force around which is described by a power law. 
Therefore, it should be possible to derive the 
same results from the complicated formulas 
representing excitation functions right above the 
threshold which have been derived from wave 
mechanics by Morse and Stueckelberg.® 





Fic. 1. The arrows indicate the directions of the electron 
impinging on and bouncing back from the molecule repre- 
sented by the circle. 


®P. M. Morse and E. C. G. Stueckelberg, Ann. d. 
Physik 9, 579 (1931). 
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A quantitative comparison of the steepness of 
this theoretical curve with the steepness of the 
experimental curves’ is meaningless. The reason 
is that the experimental curves give only relative 
intensities of spectral lines, hence relative proba- 
bilities of excitation. Thus their vertical scales 
have no absolute significance. We are only 
interested in the fact that the start of the func- 
tion from zero follows from the mechanical laws 
of impact. 


Ill. GENERAL CASE 


The principal new aspect of the present argu- 
ment is the suggestion that the same conclusion 
be applied to any collision of the first or second 
kind in which kinetic energy is used to supple- 
ment quantized energy. No quantitative state- 
ment can be expected, since the above computa- 
tion is oversimplified by disregarding transfer of 
angular momentum (assuming excitation of vi- 
bration only). Nevertheless, it is plausible to 
assume that in any such case the central collision 
will be effective with a minimum of excess 


kinetic energy of the partners. This preferred 
part played by central collisions is due to the 
fact that, at any off-center collision, a larger 
part of the kinetic energy must be preserved in 
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translation or rotation. This is another way of 
saying that the excitation function of any process, 
as well as of the simple process discussed above, 
starts from zero at the threshold and increases 
with increasing relative kinetic energy. 

(It is true that this effect of the kinetic factor 
could be counteracted by the structure factor. 
This, however, seems to necessitate very artificial 
assumptions. For example, it may happen that 
the effective impacts are confined to narrow limits 
near the direction toward the center. This restric- 
tion would eliminate the off-center collisions 
completely and permit any shape of the excita- 
tion function. We disregard assumptions of this 
type.) 

From the steep start of the excitation function 
we infer that it is unsafe to assume a uniform 
probability of transition above the threshold, as 
was implied, for example, in the measurement of 
the energy of dissociation of HO discussed in 
Section I. The practical conclusion is that non- 
equilibrium methods, although by no means 
useless, are subject to a systematic error which 
cannot be quantitatively estimated. Equilibrium 
methods, on the other hand, are free from this 
error, since a low probability of transition does 
not affect the equilibrium except by retarding the 
rate at which it is reached. 





